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Abstract 

This  paper  presents  design  considerations  and  fionda- 
mental  performance  limitations  of  phase-lock  control  loops 
which  are  used  to  coherently  combine  a specified  array  of 
lasers.  The  control  problem  is  first  investigated  in  a 
"two-laser"  control  loop  designed  to  lock  the  differential 
frequency/ phase  between  the  two  lasers  to  a specified 
reference  frequency.  An  optical  heterodyne  measurement  con- 
figuration is  used  to  determine  the  differential  frequency 
and  phase  between  the  laser  pair.  An  error  voltage  related 
to  error  between  the  desired  frequency /phase  (i.e.,  the  rf 
reference)  and  the  measured  optical  differential  frequency/ 
phase  is  filtered  and  used  to  frequency  modulate  one  of  the 
lasers  (by  electronically  changing  the  effective  cavity 
length)  in  an  attempt  to  null  the  error.  An  integro- 
differential  equation,  valid  for  the  linearized  operating 
region  of  the  loop  is  derived  in  terms  of  the  specified 
(desired)  phase  control,  the  heterodyne  measurement  noise, 
and  the  various  laser  phase  noises  (instabilities) . The 
power  spectrum  of  the  laser  phase  noise  has  two  components: 
one  due  to  quantum  limitations  (proportional  to  1/f  ) and 
the  other  due  to  random  thermal  and  acoustical  vibrations 
(proportional  to  1/f^).  The  solution  of  the  integro- 
differential  equation  results  in  an  expression  for  the  phase 

xvi 


error  variance  (the  phase  deviation  between  the  two  lasers) 

in  terms  of  the  closed-loop  equivalent  noise  bandwidth  of 

the  system,  W^.  The  variance  involves  three  terms:  the 

first  is  proportional  to  Af^/W^  where  Af^  is  the  linewidth 

of  the  combined  laser  fields  due  to  external  noise;  the 

second  is  proportional  to  Af  /W„  where  Af  is  the  combined 

q H q 

quantum  limited  laser  field  linewidth;  while  the  third,  due 
to  measurement  noise,  is  proportional  to  hfQW^  where  hfg  is 
the  energy  of  a detected  photon.  The  conflicting  effects  of 
Wjj  on  the  teinas  in  the  phase  error  variance  indicate  that 
some  optimization  of  is  desirable.  In  addition  to  noise 
effects,  the  relationship  between  the  error  voltage  filters' 
characteristics  and  and  other  system  performance  parameters 
such  as  loop  acquisition  time,  frequency  pull-in  range  and 
steady  state  error  is  examined.  A design  example  using  a CO2 
waveguide  laser  paired  with  a relatively  stable,  conventional 
CO2  laser  was  presented  and  found  to  have  an  rms  phase  devi- 
ation between  the  two  lasers  of  less  than  .1  radian  and  a 
lock  up  time  of  no  more  than  1 psec.  The  basic  pairwise  con- 
trol model  is  then  used  to  develop  control  configurations  for 
the  phase  coherent  control  of  the  laser  array.  Two  configu- 
rations are  presented:  (1)  all  lasers  are  compared  pairwise 
to  the  same  reference  laser  and  (2)  each  laser  is  locked 
, pairwise  in  sequence  across  the  array.  The  resulting  impli- 
cations of  system  complexity  and  potential  relative  phase 
errors  across  the  array  are  discussed. 
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PHASE-LOCK  CONTROL  CONSIDERATIONS 


FOR  MULTIPLE,  COHERENTLY 
COMBINED  LASERS 

I.  Introduction 

A well-known  method  of  generating  periodic,  optical 
pulse  trains  with  high  peak  powers  and  extremely  narrow 
pulses  is  the  mode-locked  laser  (Ref  1:256-266).  However, 
the  modulation  characteristics  are  largely  dependent  on  the 
laser's  cavity  parameters.  For  instance,  the  pulse  repe- 
tition rate  of  a mode-locked  laser  is  determined  primarily 
by  the  cavity  length.  Thus,  the  modulation  format  is 
inflexible. 

Conceptually,  it  is  easy  to  think  of  each  of  the  laser 
modes  as  being  generated  in  individual  laser  cavities.  So, 
a generalized  concept  of  mode  locking  could  be  the  coherent 
combination  of  the  outputs  of  an  array  of  single  mode,  but 
electronically  phase  controlled,  lasers.  Such  a configu- 
ration has  a potential  modulation  format  that  is  much  more 
flexible.  For  the  special  case  of  each  of  the  individual 
lasers  offset  from  the  next  by  a constant  frequency,  such  a 
configuration  will  result  in  the  familiar  mode-locked  wave- 
forms. More  general  phase  control  of  the  individual  lasers, 
however,  will  allow  the  synthesis  of  other  space-time  field 
waveforms . 

1 


The  key  requirement  for  such  a system  is  the  precise 
control  of  the  relative  frequencies  and  phases  of  the  indi- 
vidual lasers  of  the  array.  This  paper  presents  design  con- 
siderations and  fundamental  performance  limitations  of 
phase-lock  feedback  loops  which  are  used  to  achieve  the 
desired  coherent  laser  combination. 

The  control  problem  is  first  investigated  in  a two- 
laser  or  "pairwise"  control  loop  designed  to  lock  the 
differential  frequency/phase  between  the  two  lasers  to  a 
desired,  controllable  value.  An  integro-differential 
equation,  valid  for  the  linearized  operating  region  of  the 
loop,  is  derived  in  terms  of  the  desired  frequency /phase 
control,  the  laser  phase  instabilities,  and  other  loop  com- 
ponent noises  and  instabilities.  Expressions  for  the  power 
spectra  of  the  various  noises  and  instabilities  are  developed 
and  used  to  solve  the  integro-differential  equation  for  the 
phase  error  variance  in  terms  of  the  closed- loop  noise  equiv- 
alent bandwidth  of  the  system.  The  phase  error  variance 
expression  can  be  used  to  determine  the  optimvim,  i.e.,  mini- 
mum variance,  closed- loop  bandwidth.  In  addition  to  noise 
effects,  system  performance  parameters  (such  as  loop  acqui- 
sition time,  frequency  pull-in  range  and  steady  state  error) 
are  related  to  the  loop  filter  characteristics  and  the 
closed-loop  bandwidth  for  various  desired  control  schemes 
(such  as  step  and  ramp  changes  in  frequency) . Other  design 
considerations  such  as  loop  control  methods,  laser  field 
detectors,  laser  frequency  modulation  and  laser  quiescent 
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frequency  will  also  be  discussed.  A detailed  design  example 
for  the  two-laser  control  loop  is  presented  to  illustrate 
the  effects  of  the  various  noise  and  control  parameters  on 
the  overall  phase- lock  performance. 

Finally,  the  basic  pairwise  control  model  is  used  to 
develop  control  configurations  for  the  phase  coherent  con- 
trol of  a laser  array.  Two  configurations  are  presented: 

(1)  all  lasers  are  compared  pairwise  to  the  same  reference 
laser  and  (2)  each  laser  is  locked  pairwise  in  sequence 
across  the  array.  The  resulting  implications  of  system 
complexity  and  potential  relative  phase  errors  across  the 
array  are  discussed.  Other  array  design  considerations  (such 
as  stability  of  the  reference  laser,  laser  gain  bandwidth, 
and  beam  combination)  are  also  discussed. 
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II.  Two -Laser  Control  Loo 


The  phase  control  of  an  array  of  lasers  will  first  be 
investigated  in  "pairwise"  control  loops  designed  to  lock 
the  differential  phase  between  two  lasers  to  a desired 
value.  The  array  can  then  be  constructed  by  proper  combi- 
nation of  the  basic  pairwise  vinits.  Phase  control  is  accom- 
plished by  using  common  phaselock  techniques  which  will  be 
briefly  reviewed  first. 


Basic  Phaselock  Loop 

The  standard  device  for  estimating  the  phase  of  a 
signal  is  the  phaselock  loop  (PLL) . It  is  motivated  from 
Estimation  Theory,  where  the  maximxam  a posteriori  estimate 
of  a signal  phase  in  white,  Gaussian  noise  results  in  such  a. 
signal  processor.  Such  a model  is  valid  for  the  problem  at 
hand  (Ref  2:  Chap.  5,  p.  123).  The  basic  form  of  the  PLL  is 
shown  in  Figure  1. 


Figvure  1.  Basic  Phaselock  Loop 
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The  loop  input  reference  signal  has  the  form: 

Vj.(t)  = sin$^(t)  (1) 

where 

= (2irf^t  + + e^(t)  (2) 

The  first  term,  (2TTf^t  + f is  the  constant  portion  of  the 
reference  phase  and  represents  the  voltage-controlled 
oscillator  quiescent  frequency.  The  second  term,  ®^(t),  is 
the  modulated  term  and  consists  of  two  portions: 

\(.t)  = C(t)  + Y^(t)  (3) 

where  ?(t)  is  the  control  modulation  and  Y^(t)  is  a zero 
mean  random  process  representing  the  reference  phase  insta- 
bilities. 

When  the  PLL  is  locked,  the  phase  of  the  voltage  con- 
trolled oscillator  (VCO)  is  (ideally)  equal  to  that  of  the 
input  signal.  The  VCO  output  signal  then  has  the  form: 

Vj(t)  = cos<t^(t)  (4) 

where 

$^(t)  - (2Trf^t  + + e^(t)  (5) 

The  VCO  signal  also  consists  of  a constant,  quiescent  phase 
term,  (2iTf^t  -4-  , and  a modulated  term,  9^(t)  . 
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If  the  phase  detector  (PD)  is  assumed  to  be  linear,  its 
output,  called  the  phase  error  signal,  is  proportional  to  the 
difference  between  its  two  inputs  times  a phase  detector  gain 

v^(t)  = uK^($j.(t)  - $^(t))  (6) 

or  cancelling  the  VCO  quiescent  phase  term, 

Vg(t)  = yK^(6^(t)  - e^(t))  (7) 

where  y is  the  linear  proportionality  constant. 

The  phase  error  signal  is  filtered  by  a loop  filter, 
f (t) , which  helps  determine  the  noise  and  dynamic  perfor- 
mance of  the  loop.  The  filtered  error  signal  is 

VfCt)  = Vg(t)*f(t)  = / Vg(u)f (t-u)du  (8) 

•00 

The  frequency  of  the  VCO  is  controlled  by  the  filtered 
error  signal.  The  frequency  deviation  from  the  quiescent 
frequency  is 

0jj(t)  = 2TrKjjjVj(t)  (9) 

where  is  the  VCO  gain  in  Hz/ volt.  It  is  convenient  to 
express  these  results  in  the  frequency  domain.  Taking  the 
Laplace  transforms  of  Eqs  (7)  through  (9)  gives  respectively 

Vg(s)  - yK^(0j.(8)  - 0^(s))  (10) 
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Vf(s)  = V^(s)F(s) 


(11) 


and 


0d(s) 


Vf(s) 


(12) 


Manipulation  of  these  equations  gives  the  basic  PLL  loop 
parameter 


H(s) 


= 

0j.(s) 


KF(s) 
s + KP(s) 


(13) 


and 


l-H(s) 


0j.(s)  " ®d^s) 
Qj.(s) 


s 

s + KF(s) 


(14) 


where  K = pK  K, . The  function  H(s)  is  called  the  closed- 
m a 

loop  transfer  function  of  the  linear  loop  model.  Thus,  the 
PLL  can  equivalently  be  modeled  as  a simple  linear  filter, 
where  the  impulse  response  is  h(t)  = L~^[h(s)]  , as  shown  in 
Figure  2.  The  results  for  the  basic  PLL  will  now  be  used  to 
solve  the  two-laser  phase  lock. 


Figure  2.  Equivalent  Linear  Filter  Model  of  the 
Phaselock  Loop 
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Pairwise  Laser  Phaselock 

Ideally,  two  lasers  could  be  placed  in  phaselock  by 
letting  the  output  field  of  one  laser  be  the  PLL  input  sig- 
nal and  letting  the  other,  its  frequency  being  voltage 
adjustable,  act  as  the  VCO.  Unfortvinately,  this  scheme  is 
not  feasible  since  there  are  no  direct  detectors  for  laser 
field,  E'(r,t).  They  all  detect  the  intensity  or  power 
density  of  the  field,  |£’(r,t)|  , thus  losing  all  phase 
information.  The  extremely  high  laser  frequencies,  on  the 
order  of  10^^  Hz,  also  present  difficult  problems. 

To  circumvent  these  problems,  a common  procedure  is  to 
stun  the  output  fields  of  the  two  lasers  into  a single 
detector.  The  output  of  the  detector  will  include  a signal 
with  a phase  equal  to  the  difference  of  the  laser  phases  plus 
a low  frequency  term.  The  phase  difference  or  "beat  fre- 
quency" is  more  readily  handled  by  available  electronics. 
Typically,  one  laser  has  a very  stable  phase  so  that  the 
phase  instabilities  of  the  other  can  be  examined 
(Ref  3:509). 

Another  reason  this  procedure  is  so  attractive  is  that 
excellent  noise  performance  can  be  obtained  with  very  simple 
post-detector  processing.  This  technique  is  called  hetero- 
dyne detection,  the  results  of  which  are  examined  in 
Appendix  A. 

The  basic  PLL  of  Figure  1 can  now  be  adapted  to  lock 
the  phase  difference  between  two  lasers  to  a reference 
phase,  as  shown  in  Figure  3.  The  entire  laser  pair  and 
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Figure  3.  Two-Laser  Differential  Phase  Control  Loop 


detector  combination  acts  as  the  VCO.  The  filtered  error 
voltage  modulates  the  frequency  of  one  laser  only,  but  this 
is  translated  to  a modulation  in  the  differential  phase  seen 
at  the  detector.  The  voltage  ej^(t)  is  an  analog  control 
signal  the  uses  of  vdiich  will  be  discussed  later.  The 
important  aspect  of  the  two-laser  loop  is  that  the  differ- 
ential phase  between  the  two  lasers  is  locked  to  a reference 
phase;  the  absolute  phase  of  the  laser  pair  remains  uncon- 
trolled except  by  the  nature  of  the  lasers  themselves  and 
thus  is  free  to  fluctuate. 

VCO  Output.  It  will  be  assumed  that  the  optics  sampling 
the  laser  fields  are  aligned  such  that  the  beams  fall  collin- 
early  upon  the  detector  and  that  the  fields  are  planar  with 
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constant  amplitudes  (see  Appendix  A).  Thus  the  incident 
fields  for  the  reference  and  modulated  lasers  are  respec- 
tively: 

r 

^oCr.t)  = A^cos$q(c)  = A^Re  e (15) 

- V I" 

^^^(r.t)  = A2^cos$j^(t)  = Aj^Re  e (16) 

In  this  paper  the  Rej^  j will  be  understood  when  using  the 
complex  envelope  (exponential)  notation.  Thus  the  fields 
are  written  as 

ffo(r.t)  = A^e  ° (17) 

and 

ffj^(r,t)  = Aj^e  ^ (18) 

respectively.  The  planar  assiamption  is  reasonable  if  the 
detector  is  perpendicular  to  the  beam  at  the  beam  waist,  if 
the  detector  is  curved  to  fit  the  curvature  of  the  field,  or, 
as  is  most  likely,  if  the  detector  is  small  compared  to  the 
radius  of  curvature  of  the  field.  The  constant  amplitude 
assvimption  is  a reasonable  result  of  the  nature  of  the  laser 
instabilities  if  interest  is  primarily  in  the  field  frequency 
spectrum.  This  is  discussed  further  in  the  opening  para- 
graphs of  Appendix  B. 

The  phase  of  the  reference  laser  is 

^^(t)  =■  (2irf^t  + (j)^)  + YQ(t)  (19) 
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where  (2irf^t  + <}i^)  is  the  operating  phase  and  YqC^)  is  a 
zero  mean,  Gaussian  random  process  describing  the  measurable 
reference  laser  instabilities.  Its  statistical  description 
is  derived  in  detail  in  Appendix  B.  Since  this  laser  is  not 
modulated,  is  always  the  quiescent  frequency  fQq.  that  is. 


(20) 


The  phase  of  the  modulated  laser  is 

= (2TTf^t  + (J.p  + Yi(t)  (21) 

where  (2'irf^t  + (l)^^)  is  the  operating  phase  and  Y]^(t)  is  again 
a zero  mean,  Gaussian  random  process  representing  that  laser's 
instabilities  (Appendix  B) . More  specifically,  the  operating 
frequency  f-^,  consists  of  a constant  quiescent  frequency  ^iq 
portion  and  a portion  modulated  by  the  filtered  error  and 
control  signals.  Frequency  modulation  of  lasers  for  control 
applications  is  discussed  briefly  in  Appendix  C.  From 
Eq  (C-17) , the  instantaneous  frequency  of  the  modulated  laser 
is 

ij^(t)  = 2Tr[fj^q  + K^(Vf(t)  + ej^(t))]  + ^^^(t)  (22) 

where  is  the  modulator  gain  in  Hz  per  volt. 

The  two  laser  fields  of  Eqs  (17)  and  (18)  are  added  on 
the  detector  surface  such  that  without  loss  of  generality, 

A^  >>  Aj^.  (The  reverse  » A^  could  also  be  true  but 
would  be  included  in  the  following  with  a simple  change  of 
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subscripts.)  The  output  of  the  heterodyne  receiver  (HR)  is, 
from  Eqs  (A- 13)  and  (A- 16)  in  Appendix  A, 

Vd(t)  = V^cos$^(t)  + n(t)  (23) 

where  n(t)  is  a zero  mean  random  process  describing  the  HR 
noise.  The  phase  is 

= $3^(t)  - $Q(t)  (24) 

where 

and 

The  phase  modulation  portion  is 

e^(t)  = 2irKjj^(v£(t)  + e^(t))  + Yjj(t)  (28) 

where 

^d^^^  “ Yi(t)  - YQ(t)  (29) 

Eq  (28)  is  the  same  as  Eq  (9)  of  the  basic  PLL  results  but 
includes  a control  signal  and  an  instability  term.  The  HR 
noise  in  Eq  (23)  arises  from  the  fundamental  quantum  nature 
of  the  optical  fields  and  should  not  be  confused  with  the 
noise  arising  from  the  laser  phase  instabilities. 


2Trf^t  + + 0d(t) 


(25) 


^d  “ ^Iq  - ^oq 


(26) 


‘*’d  “ "^l  " ‘*’0 


(27) 


Using  the  quadrature  component  representation  of  the 
HR  noise,  as  outlined  in  Appendix  D,  and  assuming  that  the 
noise  spectrum  is  S3raimetric  on  the  VCO  quiescent  frequency 
fjj.  Eq  (23)  can  be  rewritten  as 

v^(t)  = V^cos$^(t)  + n^(t)cos(ai^t  + (J)^)  + (t)  sin(a)^t  + 4.^) 

(30) 

where  n (t)  and  n (t)  are  the  in-phase  and  quadrature  com- 
o s 

ponents  of  the  noise  n(t) , respectively.  Eq  (30)  is  the 
VCO  output. 

Loop  Reference  Input.  The  loop  reference  signal  remains 
the  same  as  that  of  the  basic  loop  given  in  Eqs  (1)  through 
(3) . All  derivations  will  be  made  with  respect  to  the  laser 
differential  (VCO)  quiescent  frequency,  f^. 

Phase  Detector  Output.  The  reference  and  VCO  signals, 
as  seen  in  Figure  2,  meet  in  a phase  detector  (PD).  The 
PD  will  be  modeled  as  a multiplier  with  double  frequency 
components  removed  (Ref  4:58).  The  output  of  the  PD,  called 
the  error  signal,  is  then 

where  is  the  PD  gain  and  Lp^  j denotes  the  low  pass  portion 
of  the  quantity  in  brackets.  Substituting  Eqs  (1)  and  (30) 
into  the  above  equation  results  in 

Vg(t)  - (sin<Kt)  + n*(t))  (32) 
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where 


n'(t)  = ^[n^(t)sin0j.(t)  + (t)  cos6^(t)]  (33) 

d 

and 

'l'(t)  = - <t^(t)  = 0^(t)  - 0jj(t)  (34) 

The  variable  ii;(t)  is  called  the  loop  phase  error.  As  seen 
in  Eq  (33)  the  noise  due  to  the  HR  is  independent  of  the  loop 
phase  error.  Eq  (32)  is  the  equivalent  of  Eq  (7)  in  the  basic 
PLL  results  where  y = but  is  nonlinear  in  4<(t)  and 

contains  detector  noise. 

Loop  Filter  Output.  Referring  again  to  Figure  2,  the 
error  signal  is  now  passed  through  a loop  filter  f(t)  which 
gives  a filtered  error  voltage  of 

Vf(t)  = (siniKt)  + n'(t))*f(t)  (35) 

This  is  equivalent  to  Eq  (8)  of  the  basic  PLL  results. 

Loop  Equation.  The  filtered  error  voltage  plus  a con- 
trol voltage  is  next  fed  back  into  the  frequency  modulator 
of  the  modulated  laser  changing  its  output  frequency  as 
described  in  Eq  (22) . The  instantaneous  frequency  of  the 
VCO  output  is,  from  Eqs  (25)  and  (28): 

$^(t)  = 2Tr  j^f ^ + Kjjj(Vj^(t)  + e^(t))j  + Yjj(t)  (36) 
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Substituting  into  this  expression  the  results  derived  for 
filtered  error  signal,  Eq  (35),  gives 


i^(t)  = 27rf^  + [K(sini<;(t)  + n’(t))*f(t)]  + 2TrKjjjej^(t)  + ^^(t) 

(37) 


where 


K = ttK  K,V^V , 
m d r d 


(38) 


Finally,  using  Eq  (34),  the  nonlinear  stochastic 
integro-differential  equation  for  the  loop  phase  error  is 

^'(t)  = 0^(t)  - K(sinip(t)  + n'(t))*f(t)  - 2TTKjjjej^(t)  - ;^(t) 

(39) 

This  equation  will  be  linearized  by  assuming  that  >J^(t)  is 
small  (less  than  1 radian)  so  that  sinip(t)  ==  'P(t).  Addi- 
tionally, converting  this  equation  into  the  frequency  domain 
by  taking  the  Laplace  transform  yields 

'i'(s)  = 0j.(s)  - r^(s)  - J[kF(s)(4'(s)  + N'(s))  + 2ttKjj^Ej^(s)] 

(40) 

(It  is  assxjmed  that  the  transform  of  the  sample  fxmctions 
of  the  random  processes  exist.)  This  linearized  version  of 
the  integro-differential  equation  is  graphically  depicted  in 
Figxire  4. 


Figure  4.  Linear  Two-Laser  Phaselock  Loop  Model 


Solving  for  T(s)  gives 
" [s+KF(s) 


- T— ^ 

[s+KF(s) 


2TrK^Ei(s) 


KF(s) 


s+KF(s) 


N'(s) 


(41) 


This  is  the  loop  phase  error  between  the  laser  differential 
phase  and  the  desired  reference  phase.  It  will  become  impor- 
tant when  considering  the  modulated,  steady  state  response 
of  the  system  in  a later  section. 

In  considering  an  array  of  phase  locked  lasers,  of  pri- 
mary concern  is  the  laser  differential  phase  Since 

$^(t)  = (2irf^t  + (})^)  + 0^(t),  it  is  sufficient  to  know  only 
the  modulated  phase,  9^(t).  Substituting  0^(s)  - 0^(s)  for 
4'(s)  (from  Eq  (34))  and  solving  for  ©^^(s)  gives 


®d^®^  “ s+KF(s) 


27rKj^E^(s) 


s 

r 4- 

s+KF ( s ) 

s+KF(s) 

(N’(s)  + 0^(s))  (42) 
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The  laser  differential  phase  error  can  thus  be  obtained  by 
passing  the  various  noise,  instability,  and  control  processes 
through  a set  of  linear  filters.  Eq  (42)  is  shown  graphi- 
cally in  Figure  5 for  both  Laplace  and  time  domains.  From 
Eq  (3) , the  transform  of  the  reference  modulated  phase  is 

0^(s)  = 2(s)  + r^(s)  (43) 


Figure  5.  Equivalent  Linear  Filter  Model  of  the  Two-Laser 
Phaselock  Loop  for  (a)  Laplace  and  (b)  Time  Domains 
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To  find  the  mean  and  variance  of  the  laser  differential 


phase  error,  the  statistics  of  the  noises  and  instabilities 
and  the  nature  of  the  control  modulation  and  control  signals 
must  be  known.  A first  and  second  moment  description  of  the 
noises  and  instabilities  will  be  discussed  next. 

Statistical  Description  of  Loop  Noise  Sources 

There  are  three  sources  of  noise  preventing  a steady 
value  of  laser  differential  phase.  Tney  are  the  loop  refer- 
ence phase  instabilities,  the  individual  laser  phase  insta- 
bilities, and  the  heterodyne  receiver  noise.  A second  moment 
model  for  each  will  be  presented  in  this  section. 

Loop  Reference  Instabilities.  As  defined  in  Eq  (3), 
the  phase  modulation  of  the  loop  reference  phase  consists  of 
a control  modulation  component  5(t)  and  an  instability  com- 
ponent Y^(b) . It  will  be  assumed  that  the  reference  phase 
is  very  stable  with  respect  to  the  other  noises  in  the  system, 
that  is,  ~ 0.  The  phase  modulation  term  thus  consists 

only  of  the  deterministic  control  term  so  that  the  mean  and 
spectrum  are 

E[9r(t)]  = t(t)  (44) 

and 

Sg  (f)  = S (f)  = 0 (45) 

r ' r 

(All  power  spectral  densities  in  this  paper  are  the  Fourier 
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transform  of  the  covariance  of  the  particular  stationary 
process,  not  the  autocorrelation.) 


Laser  Phase  Instabilities.  As  defined  in  Eq  (29),  the 
laser  differential  phase  instability  is  due  to  the  instabil- 
ities of  both  the  reference  and  modulated  lasers . The  mean 
and  power  spectral  density  of  the  measured  laser  instabilities 
are  derived  in  Appendix  B.  From  Eq  (B-8) , the  laser  insta- 
bilities are  zero  mean,  thus 

E[Yd(t)]  = 0 (46) 


It  will  further  be  assumed  that  the  instabilities  in 
each  laser  are  statistically  independent,  which  is  reasonable 
if  they  are  in  their  own  thermally  and  acoustically  isolated 
cavities.  Using  this  assumption,  the  spectriom  of  the  dif- 
ference is  the  sum  of  the  individual  spectra.  Thus  from 
Eq  (B-49)  is  obtained 


(f) 

^d 


+ a,) 


. ^o^co  . ^l^cl 
5 7 ' 7 7 


CO 


'cl 


(47) 


where 


and 


a 


B 


TT 


SiTlnZ 


(48) 


(49) 


The  variable  Af^  is  the  full  width  at  half  maximum  (FWHM) 
of  the  quantum  limit  of  the  laser  field  line  width  and  Af^ 
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is  the  FWHM  of  the  laser  field  line  width  due  to  all  other 
sources,  primarily  by  thermal  and  acoustic  fluctuations 
(Ref  1:304).  The  variables  f^^  and  f^j^  are  mathematically 
expedient  cutoff  frequencies  and  are  explained  in  detail  in 
the  appendix. 

Thus  the  power  spectrtim  of  the  laser  differential 

phase  noise  has  two  components:  one  due  to  quantum  limi- 

2 

tations  (proportional  to  1/f  ) and  the  other  due  to  random 
external  vibrations  (proportional  to  1/f^) . 

Heterodyne  Receiver  Noise.  The  loop  phase  noise  due 
to  the  HR  is  from  Eq  (33) 

n'(t)  = ^rn^(t)sin9^(t)  + n^  (t)  cos0^(t)]  (33) 

d 

Drawing  freely  from  the  properties  of  narrowband  noi  * listed 
in  Appendix  D,  the  mean  of  n'(t)  using  Eq  (D-2)  is  simply 

E[n'(t)]  = 0 (50) 

The  autocorrelation  of  n'(t)  is 

Rn'(t2.ti)  = E[n’(t2)n’(t3^)]  (51) 


which  when  expanded  becomes 


Rn.(t2.  V = ^ (T)E[cos(0^(t2)  - e^(tp)] 

Vd2  c 

+ R^  (T)E[sin(0^(t2)  + 

wnere  t = t2  - t^.  Since  the  noise  spectrum  of  n(t),  S^(f) , 
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is  assumed  to  be  symmetric  about  f^(Eqs  (A-14)  and  (A-17)), 

the  cross  correlation  between  n^(t)  and  n (t)  is  zero  and 

c s 

Eq  (52)  reduces  to 

R^,(t2,tj^)  = Rn  (t)e|^cosA9^(t)]  (53) 

^d 

where  ^ 

A0j.(t)  = 0^.(4)  - 

= t(t2)  - ^(t^)  + Yj.(t2)  - Yj.(tj^)  (55) 


This  result  is  strictly  not  stationary.  However,  recall 
that  it  was  previously  assumed  that  the  reference  is  very 
stable  so  that  Y^(t)  =0.  It  is  further  assumed  that  the 
desired  modulation  ?(t),  will  not  change  as  fast  as  (will 
have  a larger  correlation  time  than)  the  other  system 
processes.  So  for  time  intervals  of  interest  it  follows 
that 


A0^(t)  « 0 


(56) 


so  that 


Rn-Ctz-ti)  - ■ TT 

''d  ^ 


(57) 


or  finally 

’ A (f)  (58) 

'’d  ■= 


From  the  results  of  Appendix  A (Eqs  (A-14)  and  (A-17))  and 


Appendix  D (Eq  (D-3)),  the  power  spectrijm  of  the  loop  noise 
due  to  the  HR  is  then 


s„.(f)  = - ■ , 

‘'^oC0^d 


^1  £ 


where  q is  the  charge  of  an  electron,  n the  detector  quantum 
efficiency  in  photo-electrons  per  photon,  R is  an  effective 
HR  resistance,  is  the  active  area  of  the  detector,  and 
I/Cqi  and  are  the  magnitude  and  bandwidth  of  the  HR  BPF 
transfer  fvinction  respectively.  Eq  (59)  can  be  reduced 
further  by  using  the  value  for  stated  in  Eq  (A-15)  which 
when  substituted  gives 

Sn.(f)  ^ h ' I fl  1^01/2  (60) 

AdAin 


The  power  spectrum  of  the  loop  noise  due  to  the  HR  is  then 
proportional  to  the  energy  of  a detected  photon,  hf^. 

With  the  statistical  nature  of  the  reference  insta- 
bility, laser  instability,  and  HR  noise  now  specified,  the 
performance  of  the  laser  differential  phase  from  Eq  (42)  can 
now  be  determined. 


Steady  State,  Second  Moment  Description  of  the  Laser 


erential  Phase 


Using  the  basic  PLL  transfer  fianction  of  Eqs  (13)  and 
(14),  Eq  (42)  can  alternately  be  written  in  terms  of  the 
closed  loop  transfer  function,  H(s),  as 


0d(s)  = 


l-H(s) 


2iiK^E^(s)  + [l-H(s)]rd(s) 

+ [h(s)]  CN'(s)  + 0^(s))  (61) 


Since  the  exact  nature  of  F(s)  need  not  be  specified,  this 
equation  is  more  convenient  to  use  in  deriving  expressions 
for  the  steady  state  mean  and  variance  of  In  this 

section,  the  mean  will  only  be  derived  for  the  unmodulated 
case  (s(t)  = ej^(t)  = 0)  since  F(s)  must  be  known  if  there  is 
modulation.  The  modulated  steady  state  mean  is  discussed  in 
the  section  on  loop  filter  design.  To  simplify  the  problem, 
an  equivalent  model  for  H(s)  will  be  used. 


Noise  Equivalent  Bandwidth  Model  of  the  Closed-Loop 
Transfer  Function.  The  closed-loop  transfer  function  is 
from  Eq  (13) 

For  any  loop  filter  F(s)  considered  here,  the  following  is 
true : 

1H(0)12  = 1 (62) 

The  noise  equivalent  bandwidth  (NEB)  of  H(s)  is  defined  as 

2 

an  ideal  low  pass  filter  (LPF)  with  magnitude  (H(0)|  and 
one  sided  bandwidth  whose  area  under  the  curve  is  the 
same  as  the  area  under  lH(s)|  across  the  entire  spectrum 
(Ref  5:33).  That  is 
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7 


* • 


(63) 


iHCs)!^  ds 


or  in  Fourier  notation 


Wjj  = / lH(j2Trf)l^  df 
o 

This  is  shown  graphically  in  Figure  6. 


(64) 


Figure  6.  Noise  Equivalent  Bandwidth  Model  of  the  Closed 
Loop  Transfer  Function 


The  exact  size  of  Wjj  will  of  course  depend  on  the  exact 
nature  of  F(s).  Although  ideally  F(s)  is  only  the  loop 
filter  which  follows  the  PD  in  the  loop  (Figure  2) , this 
assvimes  that  no  filtering  occurs  elsewhere  in  the  loop.  How- 
ever, realistically,  F(s)  includes  contributions  not  only 
from  the  designed  loop  filter,  but  from  every  other  component 
in  the  loop  (e.g.,  the  response  time  of  the  laser  phase 
modulator  and  the  HR  BPF) , thus  the  component  with  the 
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narrowest  bandpass  response  will  set  the  maximum  limit  on 
the  value  of  Wjj.  In  this  model,  the  two  dominate  filtering 
soxarces  are  the  designed  filter  and  the  HR  BPF,  the  narrower 
of  which  will  determine  Wjj. 

Steady  State  Mean.  The  laser  differential  phase  error 
is  given  by  Eq  (42).  Using  the  expression  for  the  expected 
value  of  laser  phase  noise  Eq  (46);  HR  noise,  Eq  (50);  and 
loop  reference  phase  modulation,  Eq  (44),  the  mean  of  the 
laser  differential  phase  error  is 

E[0d(s)]  = 2TrKjj^E3^(s)  + Z(s)  (65) 

where 

Z(s)  = L[;(t)]  (66) 

If  the  system  is  unmodulated,  that  is  ej^(t)  = ?(t)  = 0 
then 

E[0^(t)]  = e[0^(s)J  = 0 (67) 

and  the  constant  phase  between  lasers,  from  Eqs  (2)  and  (25) 
is 

E[$d(t)]  = E[$^(t)]  = 2Trf^t  + <i)^  (68) 

That  is,  the  phase  difference  between  the  two  lasers  is 
eqiial  to  the  loop  reference  phase  when  the  system  is 
vinmodulated. 
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A similar  procedure  starting  with  Eq  (41),  results  in  a 


mean  of  the  loop  phase  error  of 


E 


s+KF(s)  1 


s+KF(s) 


ZttKjjjEj^Cs)  (69) 


Again,  if  the  system  is  unmodulated,  then 

E[i|;(t)]  = e[4'(s)]  = 0 (70) 

which  is  the  expected  result  from  Eq  (34)  since 


'l'(t)  = $j.(t)  - <fjj(t) 


(34) 


As  desired,  the  laser  differential  phase  is  matched  on 
the  average  to  the  loop  reference  phase  when  the  system  is 
unmodulated.  A measure  of  fluctuations  away  from  this  desired 
result  is  now  required.  This  will  be  found  in  terms  of  vari- 
ance from  the  mean. 


Steady  State  Variance.  Eq  (61)  can  be  slightly  modified 


as 


P ,/2ttKEt(s)  , 

0^(s)  = [l-H(s)J  + r^(s)|  + [h(s)]  (N'(s)  + 0j.(s)) 

(71) 


This  expression  shows  that  the  phase  error  is  obtained  by 
passing  the  HR  noise  and  reference  phase  modulation  through 
the  filter  H(s)  and  the  integrated  control  signal  and  laser 
instabilities  througn  the  filter  |^l-H(s)j  . Using  the  NEB 
model  of  H(s)  and  assuming  that  all  the  processes  are 
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statistically  independent  and  that  ^(t)  = ej^(t)  = 0 (no 
modulation) , the  variance  of  the  phase  error  is 


Sj^,(f)  + (f)l  df  + 2 / (f)df 

^r  - Wjj  ^d 


(72) 


The  two  terms  of  this  integral  are  illustrated  in  Figure  7 

2 

to  aid  in  showing  the  effect  of  on  the  resulting  Oq^  . 


Figure  7.  Graphical  Representation  of  Contributions  to 
the  Laser  Differential  Phase  Error  Variance  due  to  (a)  the 
Heterodyne  Receiver  Noise  and  (b)  the  Lasers'  Phase 
Instabilities 


The  total  crosshatched  area  on  the  two  graphs  is  the  vari- 
ance. In  Figure  7a,  as  previously  noted,  Boi/2  by 
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7 


definition  can  never  be  less  than  Wjj.  Thus  from  the  point 
of  view  of  the  closed-loop  transfer  function,  the  HR  noise 
n'(t)  is  nearly  white. 

Using  the  values  for  the  spectral  densities  of  Eqs  (45), 
(47),  and  (60)  and  integrating  reduces  Eq  (72)  to 


nA^A^ 


2(ao+ai) 


u: 


H 


^Cl  IT 


+ tan 


ir  , . -1 

- + tan 


(73) 


From  Appendix  B,  f^^  and  f^j^,  are  arbitrarily  small  (experi- 
mentally less  than  1 Hz  (Ref  6:184))  so  it  is  reasonable  to 
assume  that 


W. 


H 


>>  1 


(74) 


CO 


and 


W, 


H 


>>  1 


'cl 


(75) 


-1 


Using  that  result,  the  tan  terms  can  be  approximated  by 
the  first  three  terms  of  their  power  series  as 

^ cl 


tan 


•1  ^H 


TT 

T^l. 

c 


'^H 


3W, 


(76) 
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Substituting  this  approximation  back  into  Eq  (73)  finally 
yields  an  expression  for  the  laser  differential  phase  error 
variance  of 


CiWr-H 


^ A 


(77) 


where 


(78) 


and 


Cl  = 2(a^  + ap 


2{&o  + ep 
3 


(79) 


(80) 


The  variable  is  proportional  to  the  energy  of  a detected 
photon,  hf^  ; Cl  is  proportional  to  ^fq»  t^^e  sum  of  the 
quantvim  limited  laser  field  linewidths;  and  is  propor- 
tional to  Afg,  the  s\mi  of  the  laser  field  linewidths  due  to 
external  effects.  The  conflicting  effects  of  the  closed- 
loop  bandwidth  on  the  terms  in  the  phase  error  variance 
indicate  that  some  optimization  of  Wr  is  desirable. 

Two-Laser  Loop  Design  Considerations 

Several  parameters  are  important  when  considering  the 
two-laser  loop  performance.  One  measure  of  performance  is 
the  laser  differential  phase  error  variance  which,  as  was 
shown  in  Eq  (77),  depends  on  the  closed-loop  bandwidth. 

This  result,  however,  was  based  on  the  assimiption  that  the 


iT’  * / 
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loop  reference  signal  was  perfectly  stable.  Reference 
instabilities  must  ultimately  be  considered  in  the  loop 
variance  performance. 

Variance  is  not  the  only  measure  of  performance.  In 
addition,  the  loop  must  be  able  to  accurately  respond  to  the 
desired  modulation  schemes  of  the  reference  signal  thus 
placing  minimum  standards  on  the  loop  filter  design  and  the 
loop  dynamic  response.  All  of  these  design  considerations 
will  now  be  examined. 


Minimum  Phase  Variance  Bandwidth.  To  optimize  the  loop 

performance,  the  variance  of  the  laser  differential  phase 

must  be  minimized.  The  expression  for  phase  error  variance 

is  given  by  Eq  (77)  and  a plot  showing  the  effect  of  the 

three  different  terms  is  given  in  Figure  8. 

2 

To  find  Wjj  such  that  is  minimal,  take  the  derivative 
of  Eq  (77)  with  respect  to  Wjj  and  set  it  equal  to  zero.  This 
leaves 


W, 


H 


^2 


2C< 


H 


(81) 


^ 

Since  C^,  C2,  and  are  all  real,  the  only  real  solution 
to  the  cubic  equation  gives  the  minimum  variance  bandwidth  as 


(82) 


30 


Figure  8.  Closed  Loop  Bandwidth  versus  Laser  Differential 
Phase  Error  Variance 


A 2 '^3 

If  >>  C2  , which  is  typically  the  case  since,  for  most 
lasers  Af^  >>  Af^  (Ref  6:181),  then  Eq  (82)  simplifies  to 


(83) 


It  is  obvious  from  the  examination  of  Figure  7 , that  if  the 
actual  bandwidth  is  less  than  ^ the  variance  is  limited 
by  the  laser  phase  instabilities  (i.e.,  laser  phase  noise 


dominates).  If  is  greater  than  W. 


Hmin 


the  variance  is 


dominated  by  detector  noise.  Surprisingly,  it  will  be  shown 
that  in  many  conceivable  systems,  because  ^ is 
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very  large  and  thus  the  variance  is  limited  by  the  phase 
instabilities  of  the  lasers  and  is  independent  of  the  HR 
measurement  noise  (see  the  ntimerical  example  in  a later 
section) . 

I 

Loop  Filter  Requirements.  As  stated  earlier,  F(s)  in 
Figure  2 actually  includes  all  "filtering"  by  any  component 
in  the  loop.  It  will  be  assumed  that  all  of  these  compo- 
nents are  controllable  and  thus  F(s)  can  be  constructed  as 
desired. 

There  are  several  loop  filters  that  are  commonly  used 
in  phase-lock  loops  (Ref  2:23).  These  are  the  first, 
second,  and  third  order  filters  and  are  listed  in  Table  I. 

The  expressions  for  H(s)  and  Wjj  were  calculated  using 
Eqs  (13)  and  (64)  respectively. 

The  steady  state  mean  response  of  the  laser  differential 
phase  error  0|^(t)  z jt  the  modulated  case  can  now  be  examined 
for  each  of  the  different  loop  filters  under  different  con- 
trol modulation  schemes.  The  modulation  schemes  to  be 
examined  are  listed  in  Table  II  in  both  Laplace  and  time 
domains.  (Since  each  of  the  modulation  schemes  are  causal, 
they  are  unique  transform  pairs.)  The  variable  i;(t)  from 
Eq  (3)  is  the  control  phase  modulation  of  the  loop  reference 
phase  and  is  in  units  of  radians.  The  variable  ej^(t)  is  a 
control  signal  introduced  after  the  loop  filter  and  is 
measured  in  volts . 

Once  the  filter  and  modulation  scheme  are  known,  the 
mean  steady  state  response  for  the  laser  differential  phase 
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error  can  be  determined.  If  it  is  the  reference  phase  that 
is  modulated  then  the  loop  phase  error  'p(t)  becomes  impor- 
tant. From  Eq  (34) 


'l'(t)  = 9j.(t)  - 0^(t)  (34) 

and  using  Eq  (44)  (and  the  assumption  that  9^(t)  and  9^j(t) 
are  statistically  independent)  and  exchanging  terms  gives 

E[0d(t)]  = t(t)  - E[ii;(t)J  (84) 


Thus  the  mean,  steady  state  value  for  the  laser  differential 

phase  error  can  be  found  by  knowing  the  mean,  steady  state 

value  for  the  loop  phase  error  and  the  control  modulation  of 

the  reference  phase.  For  instance,  if  t(t)  is  a ramp  in 
Dt^ 

frequency,  — , and  the  loop  second  order,  then  using 
Eq  (69)  and  the  final  value  theorem,  the  mean,  steady  state 
loop  phase  is 


lim  E U(t)] 

t-*-®  ^ 


lim  se[’1'(s)] 
s-»-o 


r s"^* 

lim  s ! 

s-*-o  Ls^  + Kfe  + aK- 


i 

% 

E[iJ^(t)]ss  = D/aK  (rad) 


(85) 


(86) 


Now  utilizing  Eq  (84)  the  mean,  steady  state  laser  differ- 
ential phase  error  is 


(87) 
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That  is,  the  loop  will  follow  the  frequency  ramp  (acceler- 
ating phase)  but  there  will  be  a constant  phase  error  offset 
of  D/aK  radians.  This  procedure  is  followed  for  all  com- 
binations of  filter  orders  and  reference  phase  modulations 
schemes  and  the  results  are  summarized  in  Table  III.  Only 
the  third  order  loop  will  follow  accelerating  phase  with  no 
constant  phase  offset.  (The  first  order  loop  has  a phase 
and  frequency  offset  when  trying  to  track  accelerating  phase 
and  is  not  calculated.) 

To  calculate  the  mean  steady  state  value  for  the  laser 
differential  phase  when  the  control  signal  is  modulated, 

Eq  (65)  is  used.  For  example,  if  ej^(t)  is  a voltage  ramp, 
C^t,  and  the  loop  is  second  order,  then 

lim  [e  0j(t:)l  = lim  sE[0(s)1  = lim  s 
t-x»  ^ •*  s-*-o  ^ ■*  s-^o 

(88) 

(89) 

Thus  a constant  phase  offset  between  the  lasers  can  be 
attained  by  use  of  the  control  signal  ej^(t)  . Results  for 
other  combinations  of  loop  orders  and  control  signal  scheme 
are  listed  in  Table  III.  (The  first  order  loop  has  a fre- 
quency and  ramp  frequency  offset  when  an  accelerated  control 
signal  is  applied  and  is  not  calculated.) 

Actually,  there  appears  to  be  no  real  advantage  to 
using  ej^(t)  instead  of  ^(t)  to  control  the  laser  differential 


2ttK  C 
m r 


s^  + Ks  + aKj 
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Mean,  Steady  State  Value  of  the  Laser  Differential  Phase 
for  Different  Loop  Orders  and  Modulation  Schemes 
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This  term  has  an  additional  constant  phase  term  of  integration  dependent  on  the 
initial  conditions  of  6j(t)  (at  t = 0)  which  is  not  included  in  the  table. 


r 


phase.  The  main  use  of  the  control  signal  ej^(t)  is  in  con- 
trol of  the  modulated  laser  quiescent  frequency  f]^q*  The 
VCO  quiescent  frequency  is  from  Eq  (26) 


^d  ^Iq  ^oq 


(26) 


The  frequency  of  the  modulated  laser  from  Eq  (22)  is 

$j^(t)  = + Kjjj(v£(t)  + e3_(t))]  + ^^^(t)  (22) 

If  there  is  a dc  component  in  addition  to  any  modulated 

component  e]^^Q(j(t)  of  the  control  signal  then  Eq  (22)  becomes 


♦l(t)  = 2ir^iq  + Kjj^(Vf(t)  + e^jj^^^(t))]  + ;^(t) 


(90) 


where 


^iq  " ^Iq  ^m®ldc 


(91) 


This  will  change  the  VCO  quiescent  frequency  of  Eq  (26) 
accordingly. 

The  control  signal  could  also  be  used  to  sweep  the  VCO 
frequency  (using  a voltage  ramp)  to  "search"  for  the  loop 
reference  signal  then  release  when  the  reference  is  acquired. 
Other  dynamic  effects  of  the  loop  will  now  be  examined. 

Loop  Dynamic  Effects . The  loop  filter,  as  shown  pre- 
viously, is  important  in  determining  the  loop  performance. 
First,  the  filter  determines  the  closed-loop  transfer 
function  and  resulting  loop  bandwidth  (Table  I) . This  in 
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turn  influences  the  steady  state  variance  of  the  laser 
differential  phase  error  caused  by  the  laser  instabilities 
and  detector  noise.  Second,  the  filter  determines  the  steady 
state  response  of  the  loop  to  various  modulation  patterns 
(Table  III).  In  addition  to  these  steady  state  performances, 
the  loop  filter  also  determines  the  transient  or  dynamic 
response  of  the  loop. 

The  most  versatile  filter,  as  seen  in  Table  III,  is  the 
one  resulting  in  a third-order  loop.  The  exact  behavior  of 
a third-order  loop  is  complex  and  difficult  to  obtain 
(Ref  2:72).  However,  it  can  be  roughly  contrasted  with  the 
second-order  loop  whose  behavior  is  relatively  well  known. 

When  a third-order  loop  is  not  initially  locked  to  the  loop 
reference  phase,  it  is  not  as  stable  in  its  ability  to 
acquire  lock  as  the  second-order  loop.  The  second  order 
results  could  in  this  case  act  as  an  upper  bound  ("best  case") 
on  the -third-order  loop  performance  (Ref  7:362).  Once  locked, 
however,  the  third-order  can  out  perform  the  second-order 
loop,  most  notably  in  its  ability  to  follow  accelerating 
phase  (frequency  ramp)  with  no  steady  state  phase  error.  In 
this  case,  second-order  loop  performance  acts  as  a rough 
lower  bound  on  third-order  loop  performance.  More  specific 
results  must  be  obtained  by  computer  simulation. 

Some  important  performance  equations  for  a second-order 
loop  in  a noiseless  environment  are  listed  in  Table  IV. 

Since  they  are  valid  for  a noiseless  environment,  they  repre- 
sent the  "best  possible"  performance  in  a noisy  environment. 


39 


Equations  for  Dynamic  Performance  of  a Second-Order 
Phaselock  Loop  in  the  Absence  of  Noise 
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ACQUISITION 
SWEEP  RATE 


These  results  are  in  terms  of  radian  frequency  oi  (rad/ sec) 
instead  of  f(Hz)  since  o)  is  typically  used  in  the  literature 
when  dealing  with  dynamic  feedback  system  behavior.  The 
general  relationship  is  o)  = 2irf.  Some  results  utilize  two, 
loop-parameter  dependent  variables  called  the  natural  fre- 
quency, (0^,  and  the  damping  factor  They  are  defined  as 
follows  (Ref  7:359): 


0)^  = /iK 

(92) 

_ l./T  _ K 
" ivi  - 

(93) 

The  natural  frequency  is  the  system  phase  error  oscillation 
frequency  in  response  to  a step  in  the  input  reference  fre- 
quency when  there  is  no  damping  (5=0).  The  damping  factor 
is  a measure  of  how  quickly  the  oscillations  attenuate  as 
the  loop  acquires  lock.  Many  of  the  expressions  for  the 
same  parameter  vary  depending  on  the  reference  so  all  sources 
are  also  listed  in  the  table. 

The  two  main  areas  of  transient  behavior  are  tracking 
and  acquisition.  Tracking  parameters  assvnne  that  the  loop  is 
initially  locked  to  the  input  reference  phase.  The  hold-in 
range,  Aujj  , is  the  maximum  amount  the  reference  frequency 
can  deviate  from  the  VCO  quiescent  frequency  (HR  output 
frequency)  and  have  the  loop  remain  locked.  It  is  theoret- 
ically infinite  for  both  second  and  third-order  loops 
(Ref  7:358).  Another  important  tracking  parameter  is  the 
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frequency  step  limit,  Ao)  . It  is  the  maximum  step  in  refer- 
ence  frequency  below  which  the  loop  will  not  skip  cycles. 
Cycle  skipping  occurs  when  the  difference  between  the  refer- 
ence and  VCO  frequencies  is  large  enough  to  create  a differ- 
ential frequency  ("beat  frequency")  at  the  PD  output. 

Finally,  the  maximum  reference  frequency  ramp  that  the  loop 
will  follow  without  dropping  out  of  lock  is  called  the  maxi- 
mum  tracking  sweep  rate,  (rad/sec  ). 

Acquisition  parameters  assume  that  the  loop  is  not 
initially  locked.  The  pull-in  frequency,  AWp  , is  the  maxi- 
mum amount  the  reference  frequency  can  deviate  from  the  VCO 
quiescent  frequency  and  still  permit  the  loop  to  converge 
to  a locked  state.  This,  like  the  hold-in  range,  is  also 
theoretically  infinite.  As  the  loop  converges  on  the  refer- 
ence frequency,  the  VCO  frequency  finally  gets  close  enough 
to  the  reference  frequency  so  that  the  loop  stops  skipping 
cycles -and  settles  into  the  locked  state.  The  range  within 
which  the  loop  ceases  skipping  cycles  after  pull-in  is 
called  the  lock-in  frequency,  Aojj^.  The  associated  times  for 
the  loop  to  pull-in  and  lock-in  to  a step  change  in  reference 
frequency  are  respectively  Tp  and  Tj^.  The  pull-in  time  Tp 
includes  only  the  time  that  the  loop  cycle-skips.  The  lock- 
in  time  Tj^,  is  the  time  from  the  last  cycle-skip  until  lock 
is  acquired.  The  tracking  frequency  step  limit  Aoj  is  often 
equated  with  the  acquisition  lock-in  frequency  Aooj^  since  the 
basic  criterion  in  the  definition  of  each  is  the  maximum  fre- 
quency difference  below  which  the  loop  does  not  skip  cycles. 
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It  does  not  matter  if  this  frequency  difference  occurs  as  a 
result  of  an  input  reference  frequency  step  with  the  loop 
initially  locked  or  as  a result  of  pull-in  with  the  loop 
attempting  to  acquire  lock.  This  can  readily  be  seen  in  the 
similarity  of  the  expressions  for  Acu  and  Aw,  in  Table  IV. 

The  lock-in  time  Tj^  thus  applies  to  either  case.  Finally, 
the  maximtjm  frequency  ramp  that  the  loop  can  sweep  and  still 
acquire  lock  is  called  the  maximum  acquisition  sweep  rate, 

^max  (rad/sec^). 

It  is  important  to  note  that  these  transient  performance 
criterion  may  dominate  simple  noise  considerations.  If  for 
example,  Wjj  was  chosen  to  minimize  the  phase  error  variance 
(Eq  (82)),  it  might  not  be  large  enough  to  meet  say  the 
acquisition  sweep  rate  demanded  of  the  system. 

Loop  Reference  Instability.  The  loop  reference  phase 
^’^(t)  was  previously  assumed  to  be  perfectly  stable,  that  is, 
the  power  spectrxmi  of  the  phase  instabilities  was  zero 
(Eq  (45)).  It  is  apparent  from  careful  examination  of  Eq  (42) 
or  (61)  that  any  loop  reference  phase  instabilities  Yj.(t)  will 
be  passed  by  H(s)  and  manifested  in  the  laser  differential 
phase  error  variance. 

For  example,  if  the  power  spectral  density  of  the  refer- 
ence phase  instabilities  is  nearly  white  it  can  be  written  as 


(f) 

'r 


for  + 


where  is  a constant  and  is  a cutoff  frequency  which 

is  very  small  (Ref  8:183).  This  expression  is  substituted 

back  into  Eq  (72)  and  solved  for  phase  error  variance.  The 

additional  tern  that  arises  in  the  laser  differential  phase 

2 

error  variance  is 

2 / S (f)df  (95) 

o '' r 


which,  after  using  the  first  two  terms  of  the  series  for 
tan~^,  becomes  approximately 


N 


r 


(96) 


Not  surprisingly,  this  term  is  reduced  by  making  W„  small. 

A small  Wjj,  as  previously  discussed  can  have  an  extremely 
detrimental  effect  on  device  flexibility.  Thus  the  loop 
reference  stability  is  extremely  important  in  specific 
system  design. 

Other  possible  design  considerati "^ns  are  mentioned  in 
the  conclusion  of  this  paper.  However,  the  effect  of  any 
given  individual  parameter  (e.g.,  F(s),  Y^(t) , n' (t) , C(t), 
etc.)  on  the  overall  laser  differential  phase  0jj(t)  can 
readily  be  \ander stood  by  careful  examination  of  Figure  5 
(Eq  (42)).  The  steady  state  mean  and  variance  of  0^(t)  due 
to  unwanted  fluctuations  (detector  noise,  laser  instabilities 
and  reference  instabilities)  has  been  examined  when  the  con- 
trol inputs  were  zero.  The  steady  state  and  dynamic  changes 


in  9^(t)  with  various  control  inputs  and  filter  designs  were 
also  examined.  A design  example  for  a potentially  realizable 
two-laser  control  system  will  now  be  presented  to  illustrate 
the  effects  of  the  various  parameters  on  the  performance  of 
the  loop. 

Numerical  Design  Example 

The  specific  system  to  be  examined  here  pairs  a highly 
stable  CO2  reference  laser  with  a highly  adjustable  CO2 
waveguide  laser.  CO2  lasers  were  chosen  because  of  their 
widespread  use.  The  waveguide  laser  was  chosen  because  it 
has  a very  wide  gain  bandwidth  and  thus  can  operate  over  a 
wide  range  of  frequencies  around  the  centerline  frequency 
(Ref  9:1-33). 

The  transition  of  interest  in  this  example  is  X = 10. 6u 
which  corresponds  to: 

f^  = 2.83  X 10^^  Hz  (97) 

The  reference  laser  has  the  following  characteristics: 

Af^Q  = 6.28  X 10"^  Hz  (98) 

and 

AfgQ  = 7.77  X 10^  Hz  (99) 
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The  waveguide  laser  has  the  following  characteristics: 


^1 

* 

2.4 

watts 

(100a) 

“a 

= 

.004 

-1 

cm 

(100b) 

^1 

= 

.922 

(100c) 

R2 

1.00 

(lOOd) 

L 

= 

20 

cm 

(lOOe) 

n 

1 

(lOOf) 

and  a useable  gain  linewidth  (range  over  which  the  output 
frequency  is  adjustable  and  the  laser  power  out  does  not 

Q 

drop  below  half  its  maximum  value)  of  7.00  x 10  Hz.  From 
Eq  (B-37), 

Afcav^  = 2.85  x 10^  Hz  (101) 

It  is  further  assumed  that  the  waveguide  laser  is  operating 
far  enough  above  threshold  so  that  N2/AN  ~ 1.  Thus  from 
Eq  (B-36). 

Afqj^  = 2.00  X 10"^  Hz  (102) 


The  field  spectrum  linewidth  of  the  waveguide  laser  (includes 
both  quantxam  and  external  components)  is  measured  as 


Thus 


Afqj^  + Afgj^  = 2.12  X 10'  Hz 


(103) 


Af^j^~  2.12  X 10'  Hz 


(104) 
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As  is  true  here,  most  lasers'  linewidths  are  dominated  by 
external  disturbances  not  quantum  limitations.  With  the 
quantum  and  external  linewidths  for  both  lasers  known, 

Eqs  (48)  and  (49)  can  be  used  to  find: 

oiq  = 1.00  X 10"^  Hz  (105a) 

= 3.18  X 10’^  Hz  (105b) 

3q  = 3.47  X 10^  Hz^  (105c) 

= 2.58  X 10^^  Hz^  (105d) 

Eqs  (79)  and  (80)  are  then  used  to  obtain 

C2  = 2.06  X 10"^  sec"^  (106) 

and 

C3  = 1.72  X 10^^  sec'^  (107) 

Clearly,  Af^^^  dominates  in  the  term  C3  . The  HR  samples 
-3  2 

1 X 10  watts  (=  and  has  a quantum  efficiency  of 

n = 0.5.  Using  Eq  (78)  it  is  found  that 

- 7.51  X 10'^^  sec  (108) 

The  reference  signal  is  assumed  to  have  no  instabilities. 

If  the  HR  BPF  is  designed  to  allow  the  frequency  to  range 
over  the  entire  useable  gain  bandwidth,  then  must  be 

Q 

greater  than  7.00  x 10  Hz.  Arbitrarily,  Bq^^,  is  picked  to 
be  eqxial  to  that  amount  or : 

Bq^  = 7.00  X 10®  Hz  (109) 
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If  it  is  assumed  that  the  loop  filter  does  not  limit  the 
bandpass  of  the  loop,  then  Wjj  is  at  most  The  NEB, 

Wjj  , of  the  closed-loop  transfer  function  reasonably  must  be 
somewhat  less  than  that  and  is  arbitrarily  assumed  to  be 


Wjj  = 4.15  X 10'  Hz 


(110) 


^ ^ Ak 

Using  the  results  for  C2,  and  Wjj  the  laser  differen- 
tial phase  variance  from  Eq  (77)  is 


Oq^  = 9.99  X 10'^  rad^ 


(111) 


or  a standard  deviation  of 


Oq^  = .100  rad 


= 5.75  degrees 


(112) 


This  is  a fairly  good  result.  However,  if  is  decreased 
by  only  one  order  of  magnitude  to  4.15  x lO”  Hz  then  the 
standard  phase  deviation  will  increase  to  = 1 rad 
(57.3  degrees).  Tne  minimum  variance  bandwidth  using 
Eq  (83)  is 


Vin  = 7.7  X 10"  Hz 


(113) 


This  is  more  than  two  orders  of  magnitude  greater  than  the 
actual  Wj.  thus  the  system  is  dominated  by  laser  phase  noise. 


The  variance  at  is 


0 Imin 


8.68  X 10"7  rad^ 


(114) 
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or 

‘^Olmin  = mrad 

= .053  degrees  (115) 

which  is  the  best  noise  variance  this  particular  loop  can 
attain, 

A second-order  loop  will  be  used  in  this  example.  A 
commonly  employed  compromise  between  stability  and  speed  of 
response  for  a second-order  system  is  to  let  (Ref  2:54) 

5 “ .707  (116) 

Using  the  expression  for  Wjj  of  a second-order  loop  from 
Table  I and  Eq  (93)  to  solve  for  K gives 

K = 1.11  X 10®  sec‘^  (117) 

Again  using  Eq  (93)  gives 

a = 5.53  X 10^  sec“^  (118) 

and  Eq  (92)  results  then  in 

0)^  = 7.83  X 10^  rad/ sec  (119) 

The  waveguide  laser  is  modulated  by  a rectangular 

piezoelectric  crystal  2.5  mm  square  by  25  mm  long  with  a 

-12 

crystal  constant  of  dp^  = 225  x 10  m/volt  movinted  on  the 
rear  mirror.  Using  these  values  in  Eq  (C-14)  gives 

= 3.18  X 10^  Hz/volt  (120) 
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Using  Eq  (38)  leaves 


= 800  volts  (121) 

These  three  variables  may  be  chosen  arbitrarily  as  long  as 
their  product  is  the  above  constant.  The  detector  output 
voltage  amplitude  (Eq  (A-15))  is  most  easily  adjusted  by 
varying  the  effective  load  resistance  R.  The  loop  reference 
voltage  amplitude  is  similarly  adjusted.  The  variable 
is  the  phase  detector  gain. 

The  transient  performance  of  the  loop  is  found  by 
using  the  results  of  Table  IV.  The  tracking  frequency  step 
limit  is 

Ao)  = 2.41  X 10®  rad/sec 

S 

» 3.83  X 10^  Hz  (122) 

The  maximum  tracking  sweep  rate  is 

^max  ” 6.13  x 10^^  rad/sec^ 

= 9.76  X 10^^  sec“^  (123) 

The  acquisition  lock-in  frequency  is 

= 2.05  X 10®  rad/ sec 
= 3.26  X 10^  Hz  (124) 

As  long  as  the  changes  in  the  reference  input  frequency  are 
less  than  Ao).  (or  Ao)  ) then  the  loop  will  have  a lock-in 
time  of  no  more  than 
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= 40  nsec 


(125) 


If,  however,  the  reference  frequency  changes  by  the  full 

g 

amount  of  useable  laser  gain  linewidth,  7.00  x 10  Hz,  then 
the  pull-in  time  is 

Tp  = 722  nsec  (126) 

which  is  in  addition  to  the  lock-in  time,  giving  a maximum 
time  to  lock  of  762  nsec  for  the  system.  This  is  very  good 
and  could  allow  modulation  on  the  order  of  msec  with  the 
loop  still  in  lock.  The  maximum  acquisition  sweep  rate  is 

^max  ^ ^ rad/sec^ 

= 4.88  X 10^^  sec"^  (127) 

There  is,  however,  a steady  state  phase  error  if  this 
second-order  loop  tracks  a frequency  ramp.  If,  for  example, 
the  tracking  sweep  rate  is  at  its  maximum  allowable  value 
for  the  loop  (D^^^  = 6.13  x 10  rad/sec  ),  the  steady  state 
phase  error  from  Table  III  is  1 rad  (57.3  degrees).  A 
third-order  loop  must  be  used  to  eliminate  this  error. 

This  example  has  shown  that  even  if  one  laser  of  the 
pair  has  a large  field  spectrum  linewidth  ('*2.1  x 10^  Hz) 
the  loop  can  perform  very  well,  with  a standard  deviation 
between  the  two-laser  phases  of  less  than  6 degrees  and  a 
maximum  locking  time  of  about  762  nsec.  It  is  important  to 
note  that  the  system  performance  is  limited  by  the  laser 
phase  instabilities  (especially  that  of  the  waveguide  laser) 
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and  not  by  detector  noise.  The  key  elements  in  improving 
performance  are  to  make  Wjj  as  large  as  possible  and  to  make 
the  lasers  as  stable  as  possible. 

This  completes  the  study  of  the  two-laser  control  loop. 
The  basic  pairwise  control  model  will  now  be  used  to  develop 
control  configurations  for  the  phase  coherent  control  of  a 
laser  array. 
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III.  Laser  Arrays 


It  is  a relatively  straightforward  process  to  construct 
an  array  of  phase  locked  lasers  from  the  two-laser  control 
loop.  However,  tradeoffs  between  equipment  simplicity  and 
laser  differential  phase  performance  can  be  made.  Two  array 
systems  will  be  considered  in  terms  of  laser  differential 
phase  error  variance.  One  consists  of  each  laser  in  the 
array  paired  against  the  same  reference  laser  and  in  the 
other  each  laser  is  paired  against  the  next  in  sequence. 

For  this  discussion  the  same  quantities  that  were  used 
in  the  two-laser  problem  will  be  subscripted  with  the  appro- 
priate ntnnber  of  the  lasers  in  the  array  (n  = 0,1,2,  . . . , N)  . 
The  reference  laser  is  subscripted  "0".  Unless  otherwise 
noted,  it  is  assumed  that  the  higher  the  laser  number,  the 
higher  its  operating  frequency. 

A critical  assiomption  in  any  of  the  following  arrays 
is  that  all  of  the  loop  reference  signals  are  in  phase  with 
each  other;  that  is,  when  all  control  voltages,  e^(t) , are 
zero,  the  phase  difference  between  loop  reference  signals  is 
precisely  the  same  as  between  the  respective  laser  pairs. 
Proper  separation  (control)  of  their  phases  then  is  mandatory 
for  the  coherent  summation  of  all  laser  output  phases.  This 
requirement  will  be  discussed  later. 
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Single  Reference  Array 

A single  reference  array  is  constructed  of  N modulated 
lasers  each  paired  with  the  same  single  reference  laser,  in 
the  same  manner  as  discussed  in  the  two-laser  case,  making  a 
total  of  N lasers  in  the  array.  Thus  there  are  N separate 
loops,  one  for  each  modulated  laser  as  shown  in  Figure  9. 

The  "n^^  laser  pair"  refers  to  the  laser  pair  consisting  of 
the  reference  laser  and  laser  n.  The  "beam  sampling  optics" 
refers  to  the  optics  required  to  combine  the  laser  output 
beams  in  the  proper  requence  and  with  the  proper  alignment 
into  the  detectors.  The  "beam  combining  optics"  are  those 
optics  used  to  coherently  add  the  output  beams  and  direct  the 
sum  towards  the  target. 

It  will  be  assumed  that  the  individual  loops  are 
uncoupled;  that  is,  the  performance  of  a given  laser  pair 
does  not  depend  on  the  performance  of  any  other  laser  pair. 
The  variance  of  the  laser  differential  phase  error  for  the 
n laser  pair  is  then  given  according  to  Eq  (77)  as 


a 


^In^Hn 


/V  A 


(128) 


The  variables  62^^,  and  are  defined  by  Eqs  (78), 

(79),  (80),  and  (64)  respectively.  Since  ultimately  the 
system  adds  all  laser  beams  to  a single  output,  it  is  impor- 
tant to  know  the  variance  of  phase  differences  among  all  the 
lasers  in  the  array.  To  do  this  it  is  assiamed  that  the 
noises  and  instabilities  of  all  N lasers  are  statistically 
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independent.  The  variance  between  any  two  different  lasers 
n and  m is  then 


mn  nm  L m n J 

» E[($^(t)  - <fo(t))2]  + E[($Jt)  - <to(t))2]  (129) 

and 
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Since  a is  a measure  of  the  spread  of  the  differential 

mn 

phase  between  laser  m and  n,  then  if  n = m, 
and  thus 
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The  array  variance  is  defined  in  matrix  form  as 
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For  the  single  reference  array,  the  array  variance  is 
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If  all  the  lasers,  except  perhaps  the  reference  laser,  are 
identical  and  all  n loops  have  the  same  performance  (e.g., 
same  Wjj)  then 
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and  the  array  variance  reduces  to 
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A simplification  possible  with  this  type  of  array  is  to 
combine  all  the  heterodyne  receivers  into  one  detector  fol- 
lowed by  a parallel  array  of  bandpass  filters.  The  reference 
laser  must  be  placed  where  no  two  differential  frequencies 
are  the  same,  probably  at  the  lower  or  higher  frequency  end 
as  shown  in  Figure  10a.  If  all  the  assumptions  for  a hetero- 
dyne receiver  still  hold  for  each  laser  pair  (e.g., 

>>  , A2  , ...  , Aj^)  then  each  BPF  will  continue  to  pass 

the  same  signal  as  before  (see  Appendix  A,  Figure  A^3) . The 
main  advantage  of  this  multiple  heterodyne  receiver  is  that 
it  requires  N-1  less  detectors.  It  is  also  likely  that  this 
will  simplify  the  beam  sampling  optics  since  now  the  refer- 
ence laser  beam  is  required  at  only  one  detector,  not  N. 

There  are  some  disadvantages  however.  Where  as  before 
each  detector  needed  an  active  bandwidth  (frequency  range 
over  which  the  detector  will  respond)  only  wide  enough  to 
cover  the  frequency  difference  between  two  lasers,  now  the 
detector  must  be  able  to  respond  to  the  entire  frequency 
range  of  all  laser  outputs.  Another  disadvantage  is  that  no 
two  BPFs  can  pass  the  same  frequencies.  If  they  could,  then 
a differential  phase  of  one  laser  pair  which  drifted  or  was 
changed  too  far  would  feedback  its  signal  through  the  wrong 
loop  to  the  wrong  laser  modulator.  This  requirement  obvi- 
ously puts  a limit  on  the  flexibility  of  the  system.  Tnat 
is,  the  differential  phase  of  each  laser  pair  is  restricted 
to  a specific  range  of  frequencies  and  can  never  crossover 
into  another  laser-pair  zone  (without  changing  the  BPF  array) . 


58 


10.  Array  Design  tradeoffs  in  the  Frequency  Domain  (N=6) 


A good  compromise  to  restore  some  flexibility  and  reduced 
required  detector  active  bandwidth  might  be  to  use  two  or 
three  detectors  and  place  the  reference  laser  at  some  cen- 
tral frequency  as  shown  in  Figure  10b. 

Sequential  Array 

A sequential  array  is  constructed  of  N + 1 lasers  each 
paired  against  the  next  in  sequence,  "0"  with  "I”,  "1"  with 
"2",  "2"  with  "3",  etc.,  as  shown  in  Figure  10c,  making  a 
total  of  N control  loops . The  only  difference  between  this 
array  and  the  single  reference  is  in  the  beam  sampling  optics 
represented  in  Figure  11.  The  rest  of  the  array  is  the  same 
as  in  Figure  9 . 


Figure  11.  Representation  of  Beam  Sampling  Optics  for 
a Sequential  Array 

Again,  the  loop  reference  signals  are  assumed  to  be  in  phase 
with  each  other. 
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Assvaning  again  that  the  individual  loops  are  uncoupled, 
the  variance  of  the  laser  differential  phase  error  for  each 
pair  is  still  given  by  Eq  (77) , but  the  modulated  laser  of 
one  pair  is  the  reference  laser  of  the  next.  Thus  a change 
in  differential  phase  anywhere  will  "ripple"  down  the 
remaining  laser  pairs.  If  the  noise  and  instabilities  of 
all  N lasers  are  statistically  independent  then  the  variance 
between  any  two  lasers  is 

2 _ 2 _ 2 2 2 
^mn  “ °nn  " ‘^n(n+l)  '^(n+l)(n+2)  ^(n+2)(n+3)  ' 

•••■*■  '^(L2)(m-1)  '^(Ll)m  ' n < m (136) 

If  all  N+1  lasers  are  identical,  then 


„ 2 

°n(n+l) 

CM  r-l 

O 

O 

II 

, 0 £ n < 

N-1 

(137) 

and  the  laser 

differential  phase  variance 

• 2 
matrix  Og 

of 

the 

sequential  array  is 

0 “01 

2“oi 

3“0i 

'*“01 

0 

“01 

2“oi 

(N-l)aoi 

i 

1 

0 

“0? 

(138) 

0 

• 

• • • 

0 ! 

61 


The  further  apart  in  the  sequence  two  lasers  are,  the 
proportionally  more  is  the  variance  between  them.  This 
obviously  is  not  as  good  as  the  single  reference  array 
performance  (Eq  (135)). 

If  the  lasers  are  separated  by  a constant  frequency  Af, 
then  the  differential  frequency  for  each  laser  pair  is  the 
same.  In  this  case  a single  detector  cannot  be  used,  as  in 
the  single  reference  array,  because  each  BPF  must  pass  the 
same  differential  frequency.  However,  since  all  loops  have 
the  same  differential  frequency  Af  there  is  no  way  to 
distinguish  the  differential  frequency  of  a given  pair.  A 
possible  method  to  circumvent  this  problem  is  not  to  pair 
the  lasers  sequentially  in  frequency.  This  approach  is 
depicted  in  Figure  lOd.  With  this  configuration  the  differ- 
ential frequency  for  every  laser  pair  is  different  and  can  be 
distinguished  by  the  BPF  array.  Again,  however,  each  differ- 
ential -frequency  is  restricted  to  its  own  bandpass  (deter- 
mined by  the  BPF)  outside  of  which  it  will  modulate  the 
wrong  laser. 

A possibly  major  disadvantage  of  the  sequential  array 
is  that  if  one  laser  malfunctions,  it  will  separate  into  two 
independent,  smaller  arrays  destroying  the  device  output.  A 
laser  malfunction  in  the  single  reference  array  will  result 
in  only  one  beam  missing  from  the  device  output. 

Laser  Array  Design  Considerations 

There  are  several  important  factors  to  be  considered  in 
the  design  of  any  array.  Some  of  these  points  will  be 
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mentioned  here  but  they  are  by  no  means  a comprehensive  list 
It  may  be  helpful  to  refer  to  Figure  9 during  the  discussion 

Loop  Reference  Control.  Early  in  this  discussion  it 
was  assumed  that  the  loop  reference  signals  were  in  phase. 
Indeed  this  must  be  truf  since  the  separate  laser  phase  dif- 
ferentials cannot  be  coherent  if  the  standards  to  which  they 
are  locked  are  not. 

The  simplest  way  to  guarantee  phase  coherence  between 
reference  signals  is  to  let  them  all  be  the  output  of  a 
single  VCO.  This  of  course  places  the  restriction  that  all 
laser  differential  phases  must  be  the  same,  that  of  the  VCO. 

It  is  increasingly  more  difficult  to  control  the  dif- 
ferential frequency  and  phase  of  each  pair  separately.  This 
would  require  phase  locking  the  reference  VCOs,  or  adding  a 
control  signal  internal  to  the  loop  (see  Table  III) . Each 
of  these  procedures,  however,  would  add  to  device  complexity 

Single  Versus  Separate  Detectors.  The  advantage  of 
using  a single  detector  is  that  it  reduces  the  number  of 
detectors  required  and  thus  simplifies  the  beam  sampling 
optics.  However,  the  disadvantages  must  be  carefully  con- 
sidered. 

The  active  bandwidth  of  the  single  detector  must  range 
over  all  laser  outputs  not  just  two.  This  could  especially 
be  a problem  if  the  detector  response  is  not  uniform  (con- 
stant for  all  frequencies) . 
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Also  since  all  differential  frequencies  are  coining  from 
the  same  detector,  the  BPFs  for  the  separate  PLLs  cannot 
overlap.  Thus  each  laser  pair  is  restricted  to  a band  of 
frequencies  determined  by  the  BPF  array.  Additionally, 
limiting  the  size  of  the  BPFs  also  limits  the  bandwidth  of 
the  closed-loop  transfer  function,  W^j.  This  could  have  a 
severe  effect  on  differential  phase  error  variance  performance 
and  on  the  dynamic  performance . 

Separate  detectors,  although  more  are  required  and  beam 
sampling  optics  may  be  more  complex,  have  no  such  restrictions 
on  BPF  bandwidth,  Bq^,  and  thus  on  differential  frequency 
operation  range.  Indeed  the  BPF  may  be  a high  pass  filter 
instead,  and  cover  virtually  all  frequencies.  Also  detectors 
with  narrower  active  bandwidths  can  be  utilized  (this  will 
restrict  the  differential  frequency  operating  range,  but  not 
the  noise  performance) . 

Stability  of  the  Reference  Laser.  The  phase  locked  by 
the  basic  two-laser  control  loop  is  the  differential  phase 
between  the  two  lasers,  not  their  absolute  phases.  The 
laser  pair  absolute  phase  is  free  to  drift  throughout  the 
Laser  bandwidth.  This  is  also  true  for  the  array;  it  is  the 
’-eLatlvc  phase  among  the  lasers  that  is  locked,  not  their 
♦••oLite  phases.  If,  as  is  likely,  the  array  utilizes  fre- 
frow  as  much  of  the  laser  gain  linewidth  as  pos- 
entirely  possible  that  the  lasers  operating 
!«*•  -na  / be  forced  out  of  the  gain  linewidth  by  a 
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drifting  reference  laser.  For  those  lasers,  the  power  out 
would  then  drop  significantly  below  the  others  (or  even 
worse  cease  to  lase)  thus  having  a deleterious  effect  on  the 
array  output  field  and  phase  variance.  If  the  power  dropped 
enough,  that  particular  pair  would  drop  out  of  lock  com- 
pletely. 

Thus  the  absolute  frequency  stability  of  the  reference 
laser  is  of  paramount  importance.  There  are  many  stabili- 
zation techniques,  but  these  will  not  be  discussed  here 
(Ref  10:1015-1026).  The  key  design  requirement,  however,  is 
that  the  reference  laser  must  have  its  own  characteristics 
and  controls.  It  is  therefore  likely  that  the  reference 
laser  should  be  separate  from  the  array  lasers  and  perhaps 
not  even  used  in  the  ultimate  array  output . 

Initial  Acquisition.  The  heterodyne  receiver  output 
phase  'J’i^(t)  is  the  difference  in  phase  between  the  two  input 
lasers  (Eq  (A-12)).  However,  there  is  no  distinction  between 
which  laser  has  the  higher  frequency.  If  the  control  loop  is 
configured  to  shift  fj^  to  a higher  frequency  when  a positive 
error  voltage  is  applied  and  if  f^  is  initially  greater  than 
then  the  loop  will  drive  f^  up  to  f^  where  the  loop  will 
lock.  But  if  f^  is  initially  lower,  the  loop  will  still 
drive  f^^  higher  and  it  will  never  lock.  Therefore,  if  the 
initial  acquisition  is  to  be  assured,  the  reference  laser  must 
have  a known  frequency  relative  to  each  of  the  modulated 
lasers  at  system  turn-on.  This  is  another  good  argument  for 
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stabilizing  the  absolute  frequency  of  the  reference  laser.  If 
the  reference  laser  were  stabilized  at  one  extreme  of  the 
modulated  lasers  operating  range  even  if  a laser  pair  momen- 
tarily broke  lock,  it  would  automatically  be  in  a position  to 
relock.  The  control  voltage  could  also  be  used  to  "sweep" 

the  modulated  frequency  to  attain  lock. 

» 

Beam  Control.  The  optics  of  the  laser  array  must  fulfill 
two  functions.  First,  the  beams  must  be  coherently  combined 
and  directed  towards  the  desired  target.  And  second,  each 
field  must  be  "sampled"  for  use  in  the  control  loops.  For 
the  laser  array  fields  to  combine  with  the  desired  effect  in 
the  far  field,  several  attributes  of  the  fields  must  be 
regulated  precisely. 

Alignment  of  the  sampled  beams  into  the  detector  is 
critical  since  heterodyne  detection  requires  accurately 
parallel  and  overlapping  incident  fields.  A difference  in 
angle  between  the  two  fields  on  the  order  of  X/diameter  of 
detector  will  result  in  degraded  output  (Ref  11:186). 

Another  difficult  problem  regards  the  relative  optical 
path  lengths  of  the  beams  in  the  beam  combining  and  beam 
sampling  optics.  The  control  loop  for  each  laser  pair  locks 
the  laser  differential  phase  of  the  sampled  beam,  not  the 
actual  output  beam,  to  the  reference  phase.  So  while  the 
sampled  beams  may  be  perfectly  locked  by  the  control  array, 
the  actual  beams  directed  towards  the  target  may  not  be 
coherent  in  the  desired  manner.  To  avoid  this  problem,  for 
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each  laser  pair,  the  optical  distance  between  the  beam 
sampling  points  of  each  beam  must  precisely  equal  the  optical 
distance  between  each  sampling  point  and  the  beam  array  out- 
put combination  point  (or  be  an  integer  multiple  of  the  wave- 
length) . The  optical  precision  required  could  be  a major 
consideration  in  the  construction  of  any  array. 

These  are  only  some  of  the  array  design  considerations. 
Others  are  mentioned  in  the  suggestions  for  further  study 
following  the  summary. 
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IV.  Conclusion 


Summary 

This  paper  presented  design  considerations  and  funda- 
mental performance  limitations  of  phase  locked  feedback 
loops  which  were  used  to  coherently  combine  the  outputs  of  an 
array  of  single  mode,  but  electronically  phase  modulated, 
lasers.  The  control  problem  was  first  investigated  in  a 
"two-laser"  control  loop  designed  to  lock  the  differential 
phase  between  the  two  lasers  to  a specified  reference  phase. 
An  optical  heterodyne  receiver  configuration  was  used  to 
measure  the  differential  phase  between  the  laser  pair.  An 
error  voltage  proportional  to  the  error  between  the  desired, 
reference  phase  and  the  measured  differential  phase  was 
filtered  and  used  to  frequency  modulate  one  of  the  lasers  in 
an  attempt  to  null  the  error.  An  integro-dif ferential 
equation,  valid  for  the  linearized  operating  region  of  the 
loop  was  derived  in  terms  of  reference  phase  control,  the 
heterodyne  measurement  noise,  and  the  various  laser  phase 
instabilities.  The  solution  of  the  equation  resulted  in 
the  following  expression  for  the  laser  differential  phase 
error  variance  in  terms  of  the  closed-loop  equivalent  noise 

I 

bandwidth  of  the  system, 
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The  first  term  was  due  to  heterodyne  measurement  noise  where 
was  proportional  to  the  energy  of  a detected  photon.  The 
second  was  due  to  the  quantum  limited  noise  of  the  laser 
where  C2  was  proportional  to  the  quantum  limited  laser  field 
linewidth.  The  third  term  was  due  to  external  noise  on  the 
laser  where  was  proportional  to  the  linewidth  of  the  field 
(less  that  due  to  quantum  noise) . The  conflicting  effects  of 
Wjj  on  each  term  suggested  optimization  of  W^.  Additionally, 
the  relationship  between  the  loop  filter  and  Wjj  and  other 
system  performance  parameters  such  as  loop  acquisition  time, 
frequency  pull-in  range  and  steady  state  error  was  examined. 

A design  example  was  presented  to  illustrate  the  effects  of 
various  laser,  detector  and  desired  phase  control  parameters 
on  the  overall  performance  of  the  loop.  It  was  found  that 
the  loop  performed  very  well  for  the  specific  example  pre- 
sented with  a standard  phase  deviation  between  the  two  lasers 
of  less  than  . 1 radian  and  a lock-up  time  of  no  more  than 
1 ysec. 

The  basic  pairwise  control  model  was  then  used  to  develop 
control  configurations  for  laser  arrays.  Two  configurations 
were  reviewed:  (1)  all  lasers  locked  pairwise  to  the  same 
reference  laser  and  (2)  each  laser  locked  pairwise  in  sequence 
across  the  array.  Finally,  resulting  implications  of  system 
complexity  and  potential  sources  of  phase  errors  across  the 
array  were  discussed. 


Suggestions  for  Further  Study 

Several  of  the  assumptions  made  in  this  paper  may  not 
always  hold.  For  instance,  it  was  assumed  in  deriving  the 
differential  phase  variance  that  the  loop  reference  was  rela- 
tively stable  (Eq  (45)).  Although  treated  briefly 
(Eqs  (94)  through  (96))  a more  complete  study  is  in  order. 

It  was  also  implicit  in  the  discussion  but  never  really 
stated  that  lasers  operate  in  one  mode,  on  one  transition 
line  only.  Active  mode  competition  will  have  an  especially 
harmful  effect  on  the  far  field  but  could  also  affect  the 
loop  control.  Discrete  jumps  in  frequency  due  to  vibration/ 
rotational  level  competition  would  increase  the  phase  vari- 
ance. Both  these  effects  are  unfortunately  enhanced  by  the 
mirror  movement  which  may  be  used  to  modulate  the  frequency. 

Another  assxamption  was  that  the  lasers  in  an  array  were 
all  statistically  independent.  This  may  not  be  true  if  all 
the  lasers  are  in  the  same  structure.  The  external  portion 
of  phase  noise  in  this  case  would  be  mostly  common  to  all 
the  lasers . 

In  the  same  vein,  it  was  assumed  that  the  external 
noise  had  a Gaussian  density  (Appendix  B) . If  it  does  not, 
the  resulting  field  phase  spectrum  becomes  very  difficult  to 
determine,  but  may  be  important  if  a particular  external 
noise  source  is  persistent  (e.g.,  aircraft  vibrations  if  the 
device  is  to  be  mounted  on  an  aircraft) . 


Also  the  field  amplitude  has  instabilities  that  were 
neglected  (Eq  (B-4)).  It  will  have  an  effect  on  the  far 
field  waveform. 

The  closed-loop  transfer  function  was  modeled  as  an 
ideal  low-pass  filter  (Eq  (64)).  A more  accurate  description 
could  be  found  which  would  include  the  effects  of  other  loop 
components  (such  as  heterodyne  receiver  band  pass  filter, 
line  impedances,  modulator  crystal  capacitances  and  loop 
filter) . 

Several  assumed  linearities  are  especially  worth  inves- 
tigating. The  laser  output  power  for  instance  is  dependent 
on  where  along  the  gain  medium  bandwidth  the  laser  is  oper- 
ating. Also  the  detector  may  not  have  a imiform  response 
over  its  entire  active  range.  Further,  the  modulator  crystal 
may  not  respond  linearly  to  large  requests  for  changes  in 
frequency.  All  of  these  must  be  considered  in  the  final 
evaluation  of  performance. 

The  most  versatile  loop  filter,  the  third  order 
(Table  III),  presents  a mathematically  intractable  problem. 

A computer  simulation  of  such  a system  would  help  predict 
loop  performance  much  better  than  the  rough  estimates  using 
second  order  results  (Table  IV) . 

The  effect  on  the  array  output  of  a single  laser  ceasing 
to  lase  was  only  briefly  mentioned.  A more  thorough  study 
of  possible  array  configurations  and  their  inherent  reliabil- 
ity and  sensitivity  to  equipment  malfunctions  could  be 
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worthwhile  especially  if  the  system  is  to  be  used  in  a 
hostile  environment. 

Specific  attention  should  be  given  to  the  use  of  the 
CO2  waveguide  laser  referred  to  in  the  example.  It  combines 
a relatively  wide  useable  gain  bandwidth  with  single  fre- 
quency operation  (Ref  9:1-33). 

There  has  been  one  successful  attempt  to  control  the 
phase  of  a number  of  CO2  waveguide  lasers  (N  = 4)  and 
coherently  sxjm  their  outputs  (Ref  12:263-4).  Construction  of 
larger  arrays  with  more  flexible  modulation  capabilities 
could  now  be  considered. 
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Appendix  A 


Heterodyne  Detection 

Heterodyne  detection  is  a measurement  technique  used 
to  determine  the  phase  of  a field  relative  to  a reference 
field  while  simultaneously  suppressing  the  effects  of  thermal 
and  dark  current  noise  present  in  optical  detectors.  The 
first  and  second  moment  characteristics  of  the  heterodyne 
receiver  output  are  presented  in  this  appendix. 

The  basic  heterodyne  receiver  (HR),  depicted  in 
Figure  A-1  consists  of  mixing  (adding)  two  fields  together 
onto  the  same  detector  to  produce  a signal  with  a phase 
eqiaal  to  the  differential  phase  of  the  fields  ("beat 
frequency") . 


Figure  A-1.  Basic  Heterodyne  Receiver  Configuration 
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This  signal  is  then  passed  through  a bandpass  filter  (BPF) 
which  suppresses  the  noise  relative  to  the  signal.  If  the 
phase  of  one  field,  £’^(f,t)  called  the  local  oscillator,  is 
known,  then  the  phase  of  the  other  field  can  be  readily 
determined  from  the  detector  output. 

The  results  stated  herein  are  conditional  on  the 
following  assumptions: 

(1)  The  fields  impinging  on  the  detector  are  collinear, 
planar,  of  constant  magnitude  and  normally  incident.  They 
can  therefore  be  expressed  as 

ff^Cr.t)  = A^cos$Q(t)  = A^Rej^e  (A-1) 

and 

r 

ff^Cr.t)  = A^cos$]^(t)  = A^ReJ^e  (A-2) 

In  this  paper  the  Re^  j will  be  understood  when  using  the 
complex  envelope  notation.  Thus  the  fields  will  be  written 
as 

= A^e  (A-3) 

and 

2:3^(r,t)  = A^e  ° (A-4) 

The  phases  are 

^oCt)  = 27rf^qt  + (Dq  + (A-5) 
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and 


(A-6) 


= 2Trf^qt  + (J)^  + 9^(t) 

where  and  are  the  respective  quiescent  frequencies 
and  Q^Ct)  and  9^(t)  are  phase  modulation  terms  which  include 
both  control  and  random  phase  fluctuations. 

(2)  For  power  or  energy  calculations  only, 


f 


o 


(A-7) 


(3)  The  combined  beam  is  incident  on  area  A^  of  the 
detector,  called  the  active  area. 

(4)  The  fields,  traveling  in  free  space,  have  units  of 
volts/m//ohm  and  without  loss  of  generality. 


Ao  » 


(A-8) 


(The  reverse,  >>  A^  , could  also  be  true  but  would  be 
included  in  the  following  results  with  a simple  change  in 
subscripts . ) 

(5)  There  is  negligible  excess  noise.  (An  "excess 
noise"  term  arises  strictly  from  the  random  portion  of  the 
fields  (Ref  11:16).) 

(6)  The  BPF  is  ideal  with  impulse  response  b(t),  trans- 
fer  function  |B(j2iTf)|  , magnitude  1/Cq^  , spectral  width 
Bgj^  Hz,  and  centered  on  frequency 


f 


d 


(A-9) 
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The  transform  fxinction  is  shown  in  Figure  A-2. 


Figure  A-2.  Heterodyne  Receiver  Bandpass  Filter 
Transfer  Function 


( 


With  the  previous  assumptions,  the  second  moment  sta- 
tistics of  the  HR  output  current  ij(t)  are  (Ref  11:187) 

r 1 *inAjA^A, 

<^dM  ° -¥£  C°  ■ cos$^(t)  (A-10) 


■o^Ol 


and 


S.  (f) 


2 2 


2hf_C 


B 


01 


(A-11) 


o'^Ol 


where 


= $1  (t)  - 'I'^^Ct) 


(A- 12) 


and  h is  Planck's  constant,  q is  the  charge  of  an  electron 
' and  Ti  is  the  quantum  efficiency  of  the  detector  in  photo- 
electrons per  photon.  (All  power  spectral  densities  in  this 
paper  are  the  Fourier  transform  of  the  covariance  of  the 
particular  stationary  process,  not  the  autocorrelation.) 
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If  these  results  are  desired  in  terms  of  a voltage, 
v^(t) , multiply  the  current  i^(t)  by  an  effective  load 
resistance  R to  obtain 


E[vd(t)]  = V^cos$jj(t) 


(A- 13) 


and 


q^nAdA^ 

2“o=oi 


< Ifl  < f 


d 


+ 


(A- 14) 


where 


qnA^A^A^R 
^ ' “o‘=01 


(A-15) 


R includes  the  effective  resistance  of  the  HR  plus  any 
additional  resistance  (load)  added  to  obtain  the  desired 

Vd(t). 

These  results  are  typically  modeled  as  the  signal 

v^(t)  = E[v^(t)]  + n(t)  (A-16) 

where  n(t)  is  a zero  mean,  Gaussian  random  process  with 
spectrum  (Ref  13:1820-1823) 

Sn(f)  = (f)  (A-17) 

d 

The  HR  will  also  operate  with  more  than  one  field  added 
to  the  local  oscillator  field  as  long  as  the  resultant  beat 
frequencies  are  different  and  each  can  be  separated  by  its 
own  BPF,  The  multiple  HR  configuration  and  BPF  array  transfer 
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function  are  shown  in  Figures  A-3  and  A-4.  The  individual 
BPFs  need  not  have  the  same  magnitude.  The  output  of  each 
BPF  is  precisely  the  same  as  in  the  basic  case  and  is  given 
by  Eqs  (A- 10)  through  (A- 17) . 


Figure  A-3.  Multiple  Heterodyne  Receiver  Configuration 


Figure  A-4.  Multiple  Heterodyne  Receiver  BPF  Array 
Transfer  Function.  (Each  cell  need  not  have  the  same 
magnitude . ) 
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Appendix  B 


Laser  Phase  Instabilities 

The  objective  of  this  appendix  is  to  derive  the  power 
spectral  density  of  the  laser  phase  in  terms  of  quantities 
measurable  from  the  laser  field.  First,  the  form  of  the 
laser  output  field  and  the  nature  of  its  instabilities  will 
be  discussed.  From  there,  an  expression  for  the  spectrum 
of  the  laser  phase  in  terms  of  the  field  spectrtim  will  be 
derived.  Only  the  temporal  characteristics  of  the  field  are 
considered,  that  is  £:(r,t)  = F(t). 

The  laser  output  field  can  be  modeled  as  a sinusoid 
with  constant  amplitude  and  phase  plus  a random  noise  term 
M(t) . (This  treatment  is  taken  primarily  from  Yariv 
(Ref  1:307-318).)  The  field  is  thus  written: 

ff(t)  = Acos(2Trf^t  + (|)q)  + M(t)  (B-1) 

where  f^  and  (1)^  are  the  constant  frequency  and  phase  of  the 
field  respectively.  Fiorthermore,  the  noise  is  assumed  to  be 
narrowband  around  f^  (slowly  varying  with  respect  to  f^)  and 
can  thus  be  written  in  quadrature  (see  Appendix  D)  as 

M(t)  » c^(t)cos(2iTf^t  + (J)^)  + s^(t)sin(2T7fQt  + (fr^)  (B-2) 
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The  in-phase  and  quadrature  components,  c^(t)  and  s^(t) 
respectively,  are  Gaussian  random  processes. 

The  field  can  also  be  written  in  polar  form  as 

ff(t)  = E(t)cos(2iTf^t  + '1>Q  + Y(t))  (B-3) 

where 

E(t)  = VsnCt)  + (A  + c^(t))^  (B-4) 

and 

Y(t)  - (B-5) 


These  relationships  can  be  seen  from  the  phasor  diagram 
shown  in  Figure  B-1.  When  the  laser  is  sufficiently  above 


Figure  B-1.  Phasor  Representation  of  Laser  Instabilities 


threshold  so  that 

A >>  c^(t)  , s^(t)  (B-6) 

then  the  field  representation  reduces  to 

ff(t)  “ Aco8(2Trf^t  + <|)q  + Y(t))  (B-7) 
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where  the  phase  instability  is  now  given  by  the  approximation 

s (t) 

Y(t)  = (B-8) 

From  Eq  (B-7)  it  is  clear  that  the  spectral  characteristics 
of  the  laser  field  are  determined  almost  exclusively  by 
phase,  not  amplitude  instabilities  (Ref  6:181;  14:404).  It 

follows  from  Eq  (B-8)  that  Y(t)  is  a Gaussian  process  since 
it  is  linearly  related  to  the  Gaussian  term  s^(t) . It  will 
arbitrarily  be  given  a mean  of  zero  by  including  any  constant 
mean  terms  in  (2irf^t  + <I>q)-  However,  Y(t)  is  not  necessarily 
stationary  and  in  fact  represents  the  random  walk  of  the 
oscillation  phase  under  the  influence  of  the  various  noises 
(Ref  6:181;  14:404). 

An  assumption  about  the  stationarity  of  the  phase 
instability  Y(t)  is  now  made.  Laser  electric  fields  because 
of  their  inherently  high  frequencies  are  typically  measured 
as  a "beat"  or  relative  field  between  two  lasers  as  shown  in 
Figure  B-2  (Ref  6:182;  15:901;  16:375).  For  that  type  of 
configuration,  where  the  measurements  are  made  with  respect 
to  a standard  laser  whose  instabilities  are  negligible  com- 
pared to  those  of  the  laser  to  be  tested,  it  has  been  shown 
that  Y(t:)  is  a stationary,  Gaussian  random  process  (Ref  6:180- 
189;  15:901-902;  16:373-386;  17:277-280;  18:165-167;  19:2398- 
2407).  Such  a measiarement  technique  is  assumed  for  all  fields 
used  in  this  paper.  Thus,  any  statements  in  this  appendix 
(and  throughout  the  text)  regarding  laser  fields  or  spectra 
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will  be  interpreted  as  the  respective  measured  quantities. 
This  interpretation  is  consistent  with  the  types  of  meas- 
urements required  by  the  two-laser  phaselock  loop  discussed 
in  the  text. 


Figiore  B-2.  Typical  Arrangement  for  Measuring  the  Spectra 
of  a Laser  Field  and  Its  Instantaneous  Frequency 


Using  the  assumption  that  Y(t)  is  a zero  mean,  station- 
ary, Gaussian  random  process,  an  expression  for  the  relation 
ship  between  the  field  spectrxam  and  the  phase  instability 
spectrum  can  now  be  derived.  (This  derivation  comes  primar- 
ily from  Siegman  (Ref  6:183-184).)  The  following  cosine 
Fourier  transform  pair  will  be  utilized  (Ref  20 : 138) 

oo 

S(f)  » 2 / R(T)cos(2TrfT)dT  (B-9) 

o 

00 

R(t)  - 2 / S(f)cos(2TifT)dT  (B-10) 

o 
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where  S(f)  is  the  power  spectral  density  and  R(t)  is  the 
autocorrelation  over  the  interval  “ ^2  ~ ^1 

From  Eq  (B-7),  the  measured  laser  field  can  be  written 
in  exponential  form  as 

Y(t)) 


E(t)  = A Re 


■^j(2T:f^t  + <|>Q  + 


(B-11) 


For  notational  convenience,  the  process  x(t)  is  defined  as 


x(t)  = e 


= jY(t) 


(B-12) 


Thus  the  field  can  be  written  in  the  form 


F(t) 


x(t)e 


j(2irf^t  + <|>^)  * -j(2irf^t  + (J.^) 


+ X (t)e 


(B-13) 


where  x (t)  = e 


■jY(t) 


The  autocorrelation  of  the  field  is 


R^(t+T,t)  = E[s(t+T)ff*(t)] 


“zr  ^ 


j (2TTf  (t+T  )+({)  ) 

(x(t+T)e 

. -j(2TTfQ(t-K)-h{>^) 


+ X (t+T)e 


) (x(t)e 


j (2Trf^t-N>Q) 


+ X (t)e 


(B-14) 


Expanding  this  expression  and  neglecting  the  double  fre- 
quency terms  yields 


Rj,(t-K,t)  = -2j-|Rj^(t-l-T,t)e 


j2TrfoT  * 


+ R^(t-K,t)e 


-j2'irf„T 

n 


(B-15) 
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where 


Rjj(t  + T,t)  = E[x(t  + T)x*(t)] 


(B-16) 


( 


and 


R*(t  + T,t)  = E[x*(t  + T)x(t)]  (B-17) 


If  + T,t)  and  + T,t)  are  stationary  then  Eq  (B-15) 


becomes 


R/t)  - -5- 


r j2irf  T 


* . . -* o' 


+ R (T)e 


■j2Trf  t] 


(B-18) 


Thus  the  spectrvim  of  the  field  is 

- «c>  + <(-  « - *0>] 

This  is  a two-sided  spectrum  formed  about  f^  and  -f^.  Since 
it  is  impossible  to  distinguish  between  a positive  or  a 
negative  frequency  on  a spectrum  analyzer,  this  result  will 
be  converted  into  a one-sided  spectrum  by  flipping  the  ’ 
negative  spectrum  about  f = 0 to  overlap  the  positive  portion 
of  the  spectrvim.  If  S^(f  - f^)  is  real  and  symmetric  on  f^, 
which  is  true  if  is  real,  the  field  has  the  one-sided 

spectrvim 


S^(f)  = V 


(B-20) 


An  expression  for  S (f)  is  found  by  noting  that  the 


autocorrelation  of  x(t)  is 


R^(t  + T,t)  = E[x(t  + T)x*(t)]  * E ( 


jAY(t,T)l 


(B-21) 
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where 


AyCt.x)  = Y(t+T)  - Y(t) 


(B-22) 


Since  Y(t)  is  a zero  mean,  stationary,  Gaussian  random 
process,  the  autocorrelation  of  x(t)  can  be  shown  to  be 


(Ref  21:159-160) 


-hKAO) 

R (x)  = e 


(B-23) 


where  R^^CO)  is  the  variance  of  AY(t,x)  and  R^^Cx)  is  sta- 
tionary. Taking  the  transform  (Eq  (B-9))  gives 


S (f)  = 2 / e cos(2TTfx)dx 

^ o 


(B-24) 


Finally,  substituting  this  result  into  Eq  (B-20)  gives  the 
general  result. 


S^(f)  = A"  J e 
o 


2 ” 


cos(2TT(f-  fQ)r)dx 


(B-25) 


Next,  a relationship  between  R^^CO)  and  the  spectr\im  of 
phase  instability  must  be  found. 

To  find  R^^(O),  first  express  AY(t,x)  in  terms  of  Y(t) 


AY(t,x)  = y(t)*Y(t) 


(B-26) 


where 


y(t)  = [6(t  + x)  - 6(t)] 


(B-27) 
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The  transform  of  y(t)  is 


j2irf  T . £ 

Y(f)  = e ° - 1 = 2jeJ’^^'^sin(TrfT)  (B-28) 

so  that 

lY(f)|^  = 4 sin^(TrfT)  (B-29) 

Thus  the  following  relationship  between  the  spectra  of  Y(t) 
and  AyCt.x)  holds: 

S^^(f)  = 4 sin2(TrfT)SY(f)  (B-30) 

The  spectrum  of  Y(t)  is  not  readily  measiirable,  however  the 
spectrum  of  the  instantaneous  phase  instabilities  y (t)  = ^Y(t) 
is  (Figure  B-2) . Therefore  using  the  identity  (Ref  20:138) 


^(f) 


1 S-(f) 

7 ' 

(2Trf)^ 


(B-31) 


Eq  (B-30)  becomes 


Tsin(iTfT' 


S*.(f) 


Finally,  the  variance  of  AY^(t,T)  is 

Ri/O)  . f ^ s^^(f)df  = 2t2  /"  S;(f)df 


(B-32) 


(B-33) 


(This  result  differs  by  a factor  of  two  from  Siegman's  results 
due  to  the  definition  of  the  transforms  used  in  this  paper, 

Eqs  (B-9)  and  (B-10).)  Using  Parseval's  general  theorem 
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Eq  (B-33)  can  be  recast  into  the  form 


(B-34) 


These  are  tue  desired  results.  The  relationship 
between  the  field  spectr;am  S^(f)  and  the  phase  instability 
spectrum  S^(f)  is  given  by  Eq  (B-25)  where  R^y(O)  is  expressed 
by  Eq  (B-33)  or  (B-34)  and  (f ) by  Eq  (B-31) . These  general 
results  will  now  be  related  to  the  specific  cases  generally 
observed  in  a laser. 

The  random  phase  instability  Y(t)  of  lasers  is  composed 
of  two  components  (Ref  15:901).  One  is  an  "external"  con- 
tribution Yg(t)  due  to  acoustic  noise,  structural  vibrations, 
plasma  oscillations,  thermal  and  pressure  drifts  and  other 
environmental  disturbances.  (It  is  also  called  "technical" 
noise  (Ref  16:373),  or  "extraneous"  modulation  (Ref  17:277; 
18:165).)  The  other  element  is  a quantum  contribution 
Yq(t),  usually  much  weaker.  Quantiom  noise  is  the  quantum 
mechanical  or  statistical  limit  of  phase  fluctuations  and 
thus  sets  the  minimum  size  of  the  phase  spectrum.  This 
quantvim  noise  will  be  examined  first. 

The  generally  accepted  form  of  the  field  spectrum  due 
to  quantxim  noise  alone  for  a homogeneous  laser  operating 
above  threshold  at  atomic  line  center  is  the  Lorentzian 
lineshape  (Ref  1:154,  318;  16:374) 


.2  Af  /2tt 

Sg  (f)  - j— 

^q  ^ (Afq/2)^  + (f  - f^) 


(B-35) 
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The  variable  Af^  is  the  full  width  at  half  maximiam  (FWHM) 


of  the  lineshape.  This  linewidth  is  given  by 


1^2 


(B-36) 


where  h is  Planck's  constant,  P is  the  power  out  of  the 

laser,  N2  is  the  electron  population  of  the  upper  laser 

transition  level,  and  AN  is  the  population  difference 

between  upper  and  lower  levels.  The  variable  Af  is  the 

cav 

"cold  cavity"  bandwidth  (FWHM)  and  is  described  by 


(Ref  1:142) 


c(aj^L  - ln/I[pC^  ) 
2iTLn 


(B-37) 


where  c is  the  speed  of  light,  is  an  average  distributed 
field  loss  per  pass  (not  including  mirrors)  in  cm"^  , L is 
the  length  of  the  resonator,  n is  the  index  of  refraction  of 
the  laser  medium  and  and  R2  are  the  reflectivities  of  the 
front  and  back  mirrors  respectively. 

Assuming  that  S (f)  is  known,  S*  (f)  can  be  found 
using  Eqs  (B-25)  and  (B-33)  or  (B-34) . However  it  is  easier 
here  to  predict  a value  for  S*  (f)  and  then  show  that  it 

^q 

gives  the  desired  (f).  Therefore  the  predicted  value  for 

q 

the  spectrum  of  the  instantaneous  quantum  phase  instabilities 


S*  (f)  = 2TTAf 
fq  4 


(B-38) 


From  Eq  (B-33) 


R^^(O)  = 2ri£^t 

Substituting  this  value  into  Eq  (B-25)  gives 

00  -irAf  T 

Sg  (f)  = / e ^ cos(2ir(f-f  )T)dT 

q o 


(irAfq)^  + (2Tr(f-f^)) 


a2 

(Afq/2)2  + 


(B-39) 


(B-40) 


This  is  the  desired  field  spectrum  stated  in  Eq  (B-35). 
Thus  the  spectrum  of  the  phase  fluctuation  due  to  quantum 
noise  is,  from  Eqs  (B-38)  and  (B-31) 


(f) 


(B-41) 


The  external  noise  contribution,  as  previously  noted,  is 
due  to  a wide  variety  of  environmental  factors.  Although 
each  factor  may  have  different  statistical  properties  (i.e., 
non- Gaussian),  it  is  reasonable  to  assume  that  the  laser  field 
spectrum  resulting  from  their  combined  effects  will  be 
• Gaussian  in  shape,  and  indeed  this  is  confirmed  by  experi- 
mental evidence  (Ref  6:180-189).  The  form  of  the  field 
spectr\jmi  due  to  external  noise  only  is  thus  given  as 
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S,  (f) 


A' 

T 


V 


exp 


2ira- 


-[ 


2a  2 


(B-42) 


where  a^  is  the  standard  deviation  (rms  width)  of  the  curve 
and  is  given  by 

Og  = (8  ln2)"^Afg  (B-43) 

where  Af^  is  the  FWHM  linewidth.  Again  it  is  easier  to 
anticipate  the  value  of  S*  (f)  and  show  that  it  leads  to  the 
proper  value  of  the  field  spectrxm.  The  spectrum  of  the 
instantaneous  external  phase  instabilities  is 

4TTa2f 

S;  (f)  = -T-^  (B-44) 

^e  f^  + f"^ 

c 

where  f^  is  an  arbitrarily  small  "cutoff"  frequency  which  is 

effectively  a measure  of  how  closely  (f)  approximates  a 

e 

Gaussian  distribution.  If  it  is  perfectly  Gaussian  in  shape 
then  f^  = 0 and  (f)  becomes  an  impulse  at  f = 0.  The 

approximately  l/f2  relationship  of  Eq  (B-44)  has  been  verified 
experimentally  (Ref  6:180-189;  14:901-902;  17:373-386).  The 
corresponding  autocorrelation  is  from  Eq  (B-10) 

R.  (t)  - (2170^)^  e . T > 0 (B-45) 

This  expression  indicates  that  f^  is  difficult  to  actually 
measure  in  that  t,  the  observation  interval,  must  be  very 
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large  before  a measurable  decay  in  R*  (x)  occurs.  Since 

^e 

ideally  f is  small  R*  (x)  =a  {lua  Using  this  approxi- 

c Yg  e 

mation  for  R*  (x)  in  Eq  (B-34)  gives 
'e 


= (2TTag)2  x2 


(B-46) 


Substituting  this  value  into  Eq  (B-25)  finally  gives 
2 ” -%  (2Tra  )^x^ 

S (f)  = / e cos(27r(f-f  )x)dx 


.2 

A I?5- 

e 


(2Tr(f-fp)) 
|_  2(2Tra2)2 


21 


a2 

T 


V27ra? 


exp  - 


(f-f  ) 
' o' 


2-1 


2a' 


(B-47) 


which  is  the  desired  field  spectrum  of  Eq  (B-42) . Thus  the 
spectrum  of  the  phase  instabilities  due  to  external  noise  is, 
from  Eqs  (3-44)  and  (B-31) 


f 

^ “ -2  T Z- 

^e  at  + f ) 


(B-48) 


The  observed  field  spectrum  is  the  sum  of  the  quantum 
and  external  spectrxun  results.  Adding  the  expressions  of 
Eqs  (B-41)  and  (B-48)  gives 


1 r 2a^f 

S (f)  - S (f)  -h  S (f)  « -^'Af_  -h  ® 

* ’q  *e 


T ' T ■ n ’ 

i 'I  (f;  + f^) 

w C " 


(B-49) 
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Much  of  the  external  noise  can  be  removed  by  careful 
isolation  and  construction  of  the  laser,  but  the  quantum 
noise  is  of  course  inherent  in  all  lasers.  Typically  the 
external  noise  so  totally  overwhelms  the  quantvim  contri- 
bution that  the  shape  of  the  field  spectrum  appears,  for  all 
practical  purposes,  to  be  Gaussian  only.  It  is  only  in  the 
higher  frequencies  that  the  wider  "tails"  of  the  Lorentzian 
may  become  apparent  (Ref  16:384).  This  is  also  obvious  from 

examination  of  the  phase  noise  spectrum  Eq  (B-49) . The 

2 

qxiantum  contribution  is  proportional  to  1/f  and  the  external 
contribution  1/f^.  Thus  although  the  external  noise  clearly 
dominates  at  low  frequencies,  the  quantum  noise  persists  to 
higher  frequencies. 

This  derivation  assxmied  that  all  the  external  noises 
acted  together  to  create  a field  spectriom  that  was  Gaussian 
in  shape,  and  this  in  general  appears  to  be  true.  However 
the  Gaussian  assumption  is  destroyed  if  a single  noise 
(e.g.,  aircraft  vibrations  or  60  Hz  h;am)  is  dominant.  If 
its  specific  contribution  to  the  field  or  frequency  spectrum 
can  be  estimated,  then  Eq  (B-25)  or  (B-31)  respectively  can 
again  be  used  to  determine  its  phase  spectrvun. 


Appendix  C 


Laser  Frequency  Modulation 


Although  there  are  many  methods  of  modulating  or  sta- 
bilizing the  frequency  of  a laser  field,  the  simplest  is  by 
changing  the  optical  length  of  the  resonator  cavity  (Ref  22: 
116) . This  method  is  also  compatible  with  the  feedback 
techniques  discussed  in  this  paper.  Typically  the  change  in 
optical  length  is  affected  by  changing  the  driving  voltage 
across  a crystal  which  is  mounted  internally  or  externally 
to  the  cavity.  Both  types  will  be  briefly  discussed. 

For  a single  mode  laser,  the  axial  mode  frequencies  are 
given  by 


f 


(C-l) 


where  m is  an  integer,  c is  the  speed  of  light  and  is  the 
optical  length  of  the  cavity.  The  change  in  frequency  with 
a change  in  is 

dL, 

df  - ^ dL,  - - f (C-2) 

2LJ  ^ 


which  can  be  written  approximately  as 


(C-3) 
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The  external  modulation  technique  consists  of  mounting 
a piezoelectric  crystal  (PE)  on  the  back  of  one  of  the  cavity 
mirrors  as  shown  in  Figure  C-1. 


Figure  C-1.  External  Crystal  Frequency  Modulation 
Configuration 


Applying  voltage  v(t)  to  the  crystal  physically  changes  its 
length  according  to  the  expression  (Ref  23:115-116) 

A1  - v(t)  (C-4) 

where  1 and  Al  are  the  length  and  change  in  length  of  the 
crystal,  d is  the  width  of  the  crystal  across  which  the 
voltage  is  applied,  and  dp^  is  the  crystal  constant  in 
m/volts.  Since  expanding  the  crystal  shrinks  the  cavity, 
the  relationship  between  crystal  and  cavity  length  is 

AL  » - Al  (C-5) 

Substituting  Eqs  (C-4)  and  (C-5)  into  (C-3)  gives 

Af  - f v(t)  (C-6) 
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where 


\ = nL  (C-7) 

The  internal  modulation  case  consists  of  mounting  an 
electro-optic  crystal  (EO)  inside  the  cavity  as  shown  in 
Figure  C-2.  Applying  voltage  v(t)  to  the  crystal  changes  the 
index  of  refraction  that  the  laser  field  encounters  according 
to  the  expression  (Ref  1:347) 

Ano  = v(t)  (C-8) 


Figure  C-2.  Internal  Crystal  Frequency  Modulation 
Configuration 


where  n and  An  are  the  index  and  change  in  index  of 
o o 

refraction  of  the  crystal,  d is  the  width  of  the  crystal 
across  which  the  voltage  is  applied  and  r^Q  is  the  electro- 
optic coefficient  in  m/volt.  The  optical  length  of  the 
cavity  is 

• nL’  + n^l  (C-9) 
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where  n and  L'  are  the  index  and  length  of  the  laser  medium 
respectively.  The  change  in  optical  length  with  a change  in 
crystal  index  is 


dL,  =* 

A o 


(C-10) 


or  approximately 


AL.  = 

A o 


(C-11) 


Substituting  Eqs  (C-8)  and  (C-11)  into  Eq  (C-3)  gives 


(C-12) 


A generalized  expression  for  either  EO  or  PZT  modulation 
from  Eqs  (C-6)  and  (C-12)  is 


Af  = KjjjV(t) 


(C-13) 


where 


^m  “ ^ ^ 


(C-14) 


e » dr 


(C-15) 


for  the  external  (PE  case)  or 


’^o’^EO 


(C-16) 


for  the  internal  (EO  case) . 


The  frequency  of  the  laser  then  is  composed  of  its 
natural  or  quiescent  frequency  f^  (when  v(t)  = 0) , a modu- 
lated frequency  as  described  above,  and  any  instantaneous 
phase  instabilities,  Y(t).  Thus  the  instantaneous  frequency 
of  a modulated  laser  can  be  written  as 


Appendix  D 


Quadrature  Component  Representation 
of  Narrowband  Noise 


For  a system  operating  at  frequency  f^  with  bandwidth 
where  f^  >>  , it  is  common  to  represent  the  system 

noise  in  terms  of  quadrature  components  as 

n(t)  = n^(t)cos(2TTf^t  + (j>^)  + n^  (t)  sin(2Trf^t  + (D-1) 

where  n (t)  and  n (t)  are  the  in-phase  and  quadrature  compo- 
c s 

nents  of  the  noise  n(t)  respectively  and  <t>^  is  an  arbitrary 
phase  angle  (Ref  24:237-244).  They  have  the  following 
properties 

E[n(t)]  - E[n^(t)]  = E[n3(t)'  = 0 (D-2) 

S„  (£)  - S (f)  - L [S„(f-£j)  + S„(f+£j)]  (D-3) 

C s 

Vn  n <■>-*> 

sc  c s ^ 

where  [ ] ™eans  "the  lowpass  portion  of,"  e[  j is  the 
expected  value,  and  S(f)  is  the  power  spectral  density.  If 
S^(f)  is  symmetric  about  f^,  then  n^(t)  and  n^Ct)  are 
uncorrelated  and 


sc  c s 


(D-5) 


where  R(t)  is  the  autocorrelation.  If,  in  addition,  n(t) 

is  a Gaussian  random  process  n (t)  and  n (t)  will  be  inde- 

c s 

pendent  Gaussian  processes.  From  Eq  (D-3)  the  magnitude 

of  S (f)  (or  (f))  is  twice  that  of  S^(f) • This  is 
c s 

shown  graphically  in  Figure  D-1. 


Figure  D-1.  Spectra  of  (a)  Narrowband  Noise  and  (b)  Quad 
rature  Components  of  the  Narrowband  Noise  with  Sjj(f)  Sym- 
metric about  f^  (Ref  22:242) 
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